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ABSTRACT 


Algebraic Methods are employed for frequency response studies of 
physical systems. The transfer function of each system's mathematical 
model is composed of at least two independent variables embedded in the 
coefficients. These, together with magnitude, angular frequency, and 
phase angle are the variabie parameters to be studied. 

Computer programs are developed for graphically representing these 


variable parameters. Interpretation of these results is made to evaluate 


the utility of such methods in CEE analysis studies. 
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Te. Introduction. 

The dynamic response of linear systems can be found basically by two 
methods; experimental and analytical. Continuous advancements in tech- 
nology have precluded the feasibility of the former, relegating its utility 
to hybrid studies at best. Increased attention has therefore been focused 
on the analytical techniques available, on development of new techniques 
where necessary, or possible improvements in known techniques. A funda- 
mental dilemma to the analytic approach is the synthesis of one or more 
linearly independent differential equations which adequately describe the 
system to be studied. Whether the problem consists of multiple equations 
involving eigen values, or a single mathematical statement, the system 
limitations and characteristics are revealed by a thorough study of the 
roots of the characteristic equation. This equation is generated im- 
plicitly by matrix representation of a set of simultaneous first order 
equations, or explicitly by a rearrangement of a single higher order 
polynomial. 

Many sound advances have been made since the turn of the century in 
the realm of root location techniques for time invariant systems. As 
control engineering emerged from its embrionic stage and became embedded 
more firmly in all phases of engineering, added emphasis was placed upon 
the knowledge of systems behavior in a dynamic atmosphere consisting of 
at least one varying parameter. Among the more important achievements 
in system stability analysis are the Routh, Mikhailov, and Nyquist cri- 
teria. Frequency response characteristically has been greatly utilized as 
a fundamental system specification, and the graphical techniques of 
Bode and later Nichols have become widely used throughout the world for 


both analysis and design. Many other graphical tools have been formulated, 


including sophisticated root locus schemes involving cascade compensa- 
tion, and various forms of feedback. Inherent in all of the above, how- 
ever, has been the severe restriction of evaluating systems behavior as 
a function of only one variable parameter. The number of systems which 
truly fall within this category are in the minority, and it therefore 
becomes necessary tO employ single variable techniques on multivariable 
problems. The result is a non-optimum compromise between system sensi- 
tivity and computational reaieeiants, i.e., the calculations necessary 
to accurately predict system behavior in multivariable environments be- 
come unwieldly. 

Several noteworthy achievements have been made to partially overcome 
this difficulty. Extensive application of the analog-computer has been 
made to problems suitable for its use. The advent of the digital com- 
puter has opened a magnitude of possibilities, one of which has been the 
emergence of state space pteer.. A large potential exists in a relative- 
ly new field known as Optimal Control Theory, which is founded upon the 
maximum principal approach of Pontryagen, and Bellman's dynamic programming 
method. Stochastic processing and predicting filtering theory have become 
very influential in the field of controls due in large part to the accom- 
plishments of Wiener and Kalman. A still different approach has dealt 
with the various possible algebraic manipulations of the characteristic 
equation, hence the name ALGEBRAIC METHODS. It is this topic to which the 
following study is devoted. Basically, these methods employ techniques 
which permit engineering systems to be designed and evaluated by simul- 
taneous variation of two parameters. Extension to one or more additional 
parameter variations is possible, but again the computational requirements 


and intelligent interpretation of results become dominant factors, just as 
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in the expansion of single parameter techniques. Little imagination is 
needed to impress in one's mind the impact that this tool will have in 
optimizing two parameter design problems as its capabilities become better 
defined, and engineers achieve proficiency in its use. 

The following section presents the rudimentary mathematical theory 
as applied to electrical control systems employing two variable para- 
meters. Primary emphasis is placed upon formulating the necessary rela- 
tions for frequency response graphical representations. Subsequent sec- 
tions are devoted to the generation and evaluation of these displays from 
the standpoint of analysis and design. Frequency response and damping 
criteria are investigated by means of various forms of derivative feed- 
back. Noise limitations in physically realizable systems preclude the 
usage of derivative feedback in excess of second order in most cases, and 
for this reason the study of such stabilizing schemes is not considered 
practical. Typical examples are presented and a comprehensive evalua- 
tion of each is endeavored. Graphical representations are utilized and 
analyzed according to their respective value for information content. 

The concluding section recapitulates the most important limitations and 
advantages revealed by each graphical analysis example, and postulates 
further extensions of specific developments where applicable. A general- 
ized version of computer programs employed, with an explanatory flow- 


graph of each is included in the Appendix. 
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Live Mathematical Theory of Algebraic Methods. 

Conventional network analysis techniques rely heavily upon pole- 
zero representations in the complex frequency plane, as well as real time 
transient and frequency response curves. These methods are not conducive 
to the synthesis, or design problems, however. System parameters are 
seldom specified precisely, and a slight variation in one might correspond 
to an excessive alteration in another while still satisfying yet a third 
design criterion. The single parameter variation is seldom a rewarding 
approach to this type of problem, and it therefore becomes necessary to 
transform response curves into another reference frame where system para- 
meters are the variables. Fundamental complex variable transformation 
theory provides the necessary conversion media, based upon Cauchy's 


Principle of Argument. 


2.1 Transformation Techniques Employing Chebyshev Polynomials. 

The ensuing development is based upon three fundamental steps, and 
permits transforming the complex frequency plant into a coefficient 
plane and a parameter plane, where the corresponding variables are arbi- 
trarily designated as a> ey and alpha ( & ), beta ( 6 ) respectively. 
The first step involves choosing a contour to be transformed. The complex 
frequency damping constant zeta ( 5 ) has a radial line contour in the 
s-plane, and is selected as the transformation contour. A mapping func- 
tion must be considered next, and for this the characteristic equation is 
appropriate. The third stipulation requires a transformation variable 
which allows the mapping function to be represented by two or more inde- 
pendent equations. In this case, it was decided to separate the complex 
variable into its real and imaginary parts. Listed below are the neces- 


sary steps required to achieve this transformation. Considerable attention 
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is devoted to a complete and detailed development in order to provide 
continuity in the fundamental relationships which might be readily refer- 
enced as major divergence points for certain ensuing pursuits. The basic 
relations presented in the following diagram and adjacent equations serve 
merely to establish a consistent notation for all further developments. 
The approach, although exceedingly elementary, is considered necessary to 
avoid perpetuating existent ambiguities where similar ground rules have 


not been employed. 


5 Sem ) 





| 
5 | 
@= cos (-§) a cos == 8 — | 





2 2 
es CO = Fe wt sin OC 
i Le 2 W 
Sin aie | —cos & == es § = Wy 
= | C rm —— } a es —~ ; 9 
UU) soils. yl- 5 Wry | iOS == Wy Sip 


Substituting the relations for go and W into the top equation: 


Toe 


S=- §wptiwyyl- == wycosQ@ + J wy sin @ 


= (cos & + JsirD) 
Employing the results of the mathematician Euler, this last equation may 


ue 


be written: = Wy 


And using DeMoivre's Theorem the relation for any power of S is found 
to be: 


S* = Dy” a = Wn K(cosk@+J)sink@) (2.1-1) 
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It is now necessary to define two new variables: 
a \ 
di. (-$) = oc kk B= cos (k Cag (ale) 


" ae oe we Weres: F's) 
lat §) sin @ Sin laescoeye ~' (= §) 


Substitution of these equivalences into (2.1-1) yields: 


fe, a ce 


Since engineering problems involving stable systems are concerned with 
positive values of zeta, it is advantageous to employ the following 


relations which are obtained by trigonometric identities: 


eh cee an i veale’,). cae al alent) 


All characteristic equations can be reduced to a general form which is 


herein defined as: 


2 n=! nN 
Poa, et, S a5 S 4 ore Garey ay, S = (- 


Rewriting this in summation notation, and employing the identities 


developed above, we can write: 


F(s) = 


0 


IMs 
oe 
\{ 
aM 
S 
9 
x 
“e 
Nx 
BX) 
= 
bs 
{I 
© 


nN k | _ 7 Kk 
po) = SO ml pe Cos hile tel Somoup tug eSinekcO =a 
kK=0 


\ 
K=O 
n vn 
ae K : Fs k Rack 
Fed = 2 ay wy CNT, CS) td a, Wy Cl) Uy (5) =O 
A polynomial may be set to zero when its orthogonal components are 


independently zero; hence: 


K 


z Boel. ei ‘la? (¢) =o) (2.1-2) 


a, Wy (1)! Uy (E) =O (2.1-3) 


a Ms 


0 


It would now seem plausible to inspect several values of these new vari- 
ables Tag (¢) and U yn ¢) , in order that any discernable character- 


istic traits might be utilized. 


,¢5)= (-1)" 7-9) = Te (-§) = cos(08) = eos(o)=! 


—2e-—CM TC |, Ce) = cost(cos woe & 


a SSintee)... ie = 
sin (6) Sin (6) 


ies) — -)atee: 


2 . sal 
a _¢) a sinleos (-§)) 
U, eo) = ( ) U, ( a Sia (Gos | Cop) = 


Comparison of these polynomials with a table of Chebyshev functions will 
reveal their identity. Accordingly, the following recurrence relations 


may be employed for complete tabulation. 
TecennaleS ae — 2 Seale (6 teen Get ee =a 


Soe) =—2. < Vets) +c ene’) =— O 


It can be shown that [,({) and U,({) are related by the following: 
Tilt i isee Selly Calder TA pe (GD 


Substituting this into (2.1-2) yields: 


Fey Wy Cl) "(Uy i 6)] — Z ay wyKCl)"[U,()]I=O 2-H 


which can be rewritten as follows: 


i: 


~ y . 
$2 ay wy (NE Uy QI— 2 ay wy“ (-1)" Cup (SIE 0 


Factoring a (-1) out of (2.1-3) and multiplying that result through by 


( § ) in a separate step yields the following sequential equations, so 


arranged for easy inspection. 


:. ay Wy," (~1)" U-(§) =O= sls. oa Wy" (-1)" Uy (8) | 


This last form of (2.1-3) is identically equal to the first term of 


(2.1-4). Therefore the two fundamental relations become: 


M3 


(~1) " ay Wy ea) (2.1-5) 


NN 


=0 


7 M2 


en" Qa Ww, Uy CoO (2.1-6) 
The utility of (2.1-5) and (2.1-6) is not obvious, but their judi- 
cious use results in a powerful tool. Employment of the bilateral trans- 
formation allows the natural angular frequency ( W), ) to range from 

minus infinity to plus infinity, thereby eliminating it as an unknown 
variable per se. Choosing a specific value for zeta ( g ) allows the 
values of U (¢) to be obtained by the preceding method or by refer- 
ence to precomputed tables [1]. Thus any two values of Q&y can be al- 
lowed to vary, thereby formulating two independent equations in terms of 
two variable coefficients of the characteristic equation. Mitrovic [2], 
who pioneered the above development, chose the two lowest order coef- 
ficients, Q, and Q, , for his study. Elliott, Heseltine, and Thaler, [3] 
and summarily Siljak [4] extended Mitrovic's approach to additional co- 
efficients. A brief summary of the coefficient plane derivation appears 


on the next page. 
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2.2 Coefficient Plane 

For the generalized form of the characteristic equation, let two 
variable coefficients be defined as Q, and Q y: The order of the 
equation ( 7 ) is no smaller than xX , and X is greater than y. 
The value of y is non-negative. (mn =X > Y>0) Transposing the @, 


and ay terms in (2.1-5) and (2.1-6) yields the following results: 


Y 


fy a, Ww, Hy. (6) + CN) aw, Uy, Cs) 


Y K K 
> Be Fee may 
KSo 
K#X 
KAY 


he ae) 


yk ow, We) + Cae, Ua 


= z (-1) "ag Wy" U C$) 


kK #Xx 
kfy 


Solving the above equations simultaneously for ay and Ay yields: 


aut oc. Uy UO) =U y CE) ee ee 
k=0 KU, (8) U, Bey) Ug (i) eee 
k Hx 
k#y 
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nN =, aN “Tet SN 
are ay 4 a apne) U, (¢) Vy (SU, alee me, 25 
UL, Uy, (9) —4,., (DU, 5) 


It will now be shown that: 
i ears) oc Hove eee (5) —U,., Co) Ce) 2. 2ew 


Let (a 


Uy (9) = U,(4) US) — UZ) Uy, (9) 
From the preceding Chebyshev recursion relations: 

U(§) =| 

Cre en 6) eel) — Chava 


Therefore: 


bh, CS inns ) eS laeenS) (2.2-4) 


| 


But it has been shown that: 


vee CS) =e : ey) (2.2-5) 
Substituting (2.2-5) into (2.2-4), we see the general recursion formula 
as previously stated for Chebyshev polynomials: 
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Ws a NE ale’) «Manes Gale @ 


It can be shown by inductive reasoning that the above relation exists 


for all values of kK [5]. Therefore we can also write: 


May, C8) = “Uys SU ol 8). ally” 65) Ua, 68) 


Thus the general solution for two variable coefficients becomes: 


n K-X (=X 

— oe (—}) wy Cm Un-y (8) _ 
— Ux-y ( §) 
kK séx 


KA y 


7 K-Y —x-y — 
hy OE eR Pt, Bren (8) 
= = & S) 
Ax 
Hy 


XK WX 


A graphical display of the characteristic equation is often formed with 


CL x as one axis and Ay as the other axis. Such a plot is termed a 


coefficient plane. With sufficient data and interpolation, it is pos- 


sible to derive a very satisfactory design in this manner for systems whose 
characteristic equations possess the desired variables in coefficient form 
An excellent treatise on the applied coefficient plane can be found in 


the NASA Contractor Report by Thaler, Siljak, and Dorf [6]. 
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2.3 Parameter Plane, Linear Case. 

The coefficient plane provides a means of working with two variable 
coefficients of a polynomial, as long as each of these variables does 
not appear in more than one coefficient. It can be recalled that each 
coefficient within a polynomial is merely the sum of the product of the 
roots taken 7 at a time, where Y equals the order of the equation 
minus the order of the particular coefficient's independent variable. 
Therefore severe restrictions are encountered when attempting to use 
this method for simultaneous, but independent variation of any two roots 
of a polynomial. Thus the coefficient plane method does not include 
the more general class of problems which have the desired variables em- 
bedded within several coefficients. Siljak [7] postulated that this 
class of problems could be handled in a manner similar to that of the 
coefficient plane, provided these variables appeared linearly. Desig- 
nating the variable parameters as alpha ( c ) and beta ( G ), one 
can rewrite any coefficient of a general polynomial as follows: 

Aavy—= b,A+C, 8 + dy 
Using this expanded form, the two independent equations (2.1-5) and 


(2.1-6) become: 


a 
POD 
nis 


a 
| 
a” 

aS 
=~ 
c~ 
> 
Sn, 
uw 
., 
4 
| 
oS 
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Kk 
n “x« 


n K . 
a= Cl) & w ( 
k=6 ; 


Y) K 
Sin) eee / 
a ) O 
“ Ww <= 
2&1) K n ue 
This rather unwieldy notation will be slightly modified by the follow- 


ing substitutions: 


Let: 


elaine ll parietal! hale 


ae ee (-1)* d , Wh,” Uy, (§) 


Kero 
Employment of these substitutions reduces (2.1-5) and (2.1-6) to the 


following form: 
BP CS Wied. tua Chal Suge ) ata) (5, Olle enmma 


x ma |) Wy, } = Cue we = Ds (S$, uw.) =O 


21 


These two equations are quite easily solved for aipha and beta by Cramer's 
rule, or by employing the algebraic method of elimination. The resultant 


solutions are: 


C, C8, Wn) Do ($n) — Co (§) wn) 08, %) 





Se Ge j : (2.3-1) 
feo ( 50g) C eee wy uc, (eee 
mee? > Se Wy) By ($,"n) =~ Do ( Sy Ww), ) B C$, nr) 

6 en a ne UR) 
B, C$, Wy) C, (£ , Wy) om By (¢)n) Cen 


The stability variables most commonly employed by control engineers 
are the damping ratio, zeta, and the undamped natural frequency, (Wy). 
Most of the current literature places great importance upon these terms, 
and interprets system performance by correlation with the standardized 
second order system. In some cases, however, the rectangular complex 
frequency plane configuration consisting of the imaginary frequency 
sigma ( 0” ) and the real frequency omega ( LL’ ) is more useful than 
the polar coordinate system with damping ratio, zeta ( ¢ ) and the 
undamped natural frequency ( Wy ). In second order system parlance, 
sigma represents the rate of exponential decay of the response signal, 
while omega is the angular frequency at which this signal oscillates. 
The necessary steps for representing (2.3-1) and (2.3-2) in terms of 
sigma and omega involve the use of the binomial expansion. A brief 
outline of this method is now presented. Recalling the basic relations 
between 0 , W , € and Wy in 2.1, it is possible to rewrite (2.1-1) 


as follows: 
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a = . He a 
Se alg We = Co ei a Xe et jee 
where the individual components of Xx and yY, are seen to have a 


particular recurrence pattern as shown by the following tabulation: 


7 

5) =o +jw 

ge 2 UL 

See = Bo ee) Soe — oe) 


G9 ob oetwt ew +) (4 oOo tw -4¥erw?) 


Go we - 10 Oi ae 5S aed” da So Ww Oo w ee 


From which: 


Me | 

—> oO 

a o 2 —-w 

X2 mo 2 — 30 w* 

ot ~ 6 we wit 
fo == © 

mee 

7 20 Ww 

¥2 = Regen tengyy ~w 

_ = 4r3zw -4orw? 


which can be consolidated into binomial expansion form for the general- 


ized case. 


g Kae] (she" = ge ee paca oe 4 
Pei) oe salon ae fa --- | = £— 
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s k 2 7[ 0" Ee aa! pee 7 
ss j lg ulin nutes ma ioe 4 Gelmeak© 


.) x | . 
where | x = “Ce-HybOx! fOr fees X 
and ( : ) —= whenever {x < X 
Recalling the generalized form of the characteristic equation F)=£4, $ 


Kz0 
we can again separate § a into its perpendicular components XK and VY, 


such that z Gay K XK represents the summation of the real terms, 
K=o 
er 4x Ya denotes the summation of all imaginary terms. Thus: 


TNS 
ta 
| 

M3 

Wx 


in aor fom 
(Ve wae ome ul 


Ce 

O K=0 j=0 wl 

yn a ee ee k K-2t+/ 2+! 
z Q ic YK = 2 ze ay [ ic I) Coins) a tad J 


In order to effectively utilize these equations for computer applica- 


tions, it is necessary that the recursion relation be explicitly formu- 


lated. This can be achieved by observing the following: 
¢ 


multiplying these two equations; 
(S-o-jJW)(S$-o7+j Ww) =O 
Expanding this, we have: 


Z a 
ely ee ai 


Pre-multiplying this by ¢ we may write 


5 ‘ [ G2 = 2 5 eee fe Uae = © 
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which, upon expansion, becomes 


ae 
oO See pee See a = OO 


This may be rewritten in terms of orthogonal components where 


Vr Die i 
5 — Oe a ig! ie 
K+) __ : 
S — eee oe "KE 
5s = Ning an ve 


Aieyep + J Pues PaO pees, tI Y%pe)) 
# (a+ w*)(X, + J yoo 
Separating this into its real and imaginary parts results in 
2 a 
Mein — 2 OOK spaieye E OES 


Oe 2a _ 
7ks2 ~~ 7 aa, oe Ce aay ee 
Noting that 0 is xX, , andwW is y, , we can make these substitutions, 


resulting in the following explicit recurrence relations: 
mn ~2X%,%x +(X 27+ y7) Xy =O 
K+ ike! f n 
' — ae \ 5 Vv _ 
Vise me— 2 Kn, Nm SOK. 7) Fe a 


Thus one can solve a given coefficient equation with coefficients of the 
form Q, = yeah +*G, 0 +d. for the parameters alpha and beta in terms 
of the variables sigma and omega. The method is the same as for the case 


involving the variables zeta and w,, . The two independent equations 


become: 
n 2 < 
oa i oe 
Oe bik + 2 6m, + 2 A, ee =O 
K=0 K=0 Nzo0 
n n Nn 
Oe =a; Yeer fF 2 Ge "K) & ce Ag yr = oO 


which may be written in abbreviated form as 


Ob, Coo) - PC, ae) + lo, 0) sme oO 


0 By (Giga )eettieeC suladagitndd. Wr D,(o,w) = O 


where: 
uy n 

eo Za Ban eo a é be Ny 
# n 

a oa Cx Xe a ee CK Nx 
/ PAs 

D= 2 dk Xx Do= 2, 4k Me 


As in the preceding development, the solutions for alpha and beta take 


the form: 

= Cito, w) Dy low) — Cz low) D, (jw) | (24233) 
BiG, w) C2 (au) — Br (aw) C, Cg w) 

e= Duley Ww) Bo. (a; w) — Pao) Bi (ow), @: 3-4) 


B, (ow) Cy (Ow) — Poa (75 w) €, (Gx) 


These equations are seen to be of the same form as (2.3-1) and (2.3-2). 
Thus the method of solution is identical to the one using Chebyshev 
se tvmend Seen with only the recurrence relation and the solution variables 
being different. 

An analogous approach to this problem does not involve the Chebyshev 
polynomials method at all, as can be seen by the following procedure 


which starts by employing (2.1-1) directly: 


N 
K 
F(s)= 2 a, S$ = 22, Wy © cos ké@ 
kK=0 K=O 
Som K 
oe) a Ak Wy, Stn ke = O 
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Separating this polynomial into its orthogonal components we may write: 


ery = Oe om cos K@Q =O 


| 


ImfFO)}= 2 ay wy" sin kO =O 


K=O 


Now substituting the integral components of A, into each of the above 


yields the following general terms: 


n nN 
i k 
i Piet CCS KOT Be Gi ctu ecotulge 
K=6 K=0 
n iL 
pas i oe Ft ae 
K=0 
nN n K 
a 2 ry Ww, sin fo + 8 Ss Cau gin k@ 
a ‘ ( ey) 


Solving these equations simultaneously for alpha and beta, we have: 














Z k 
2  Caioe coskOd dy Wy, ne “£6 Cy Wi, eee coskO 
x = -—< ~----———- —-—- eo a = ji (2.3-5) 
Ss by Ww, i kO Zc, Wy, “sinkO —Z by Wy, “sin KO BC yu sink® 
keo K=o0 A= 
k 
sub Wy, kK oskO& b Wy ‘sink@ -2 dw, i al Wy, cask 
B= a a (2.3-6) 
S be Wy,  cosk® 26, WwW, * sink@- Fale aus sink Z Cy Wy Koos KO 
kK=C =0 
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These two alternative parametric solutions are seen to be no less compli- 
cated in form than their Chebyshev polynomial or binomial expansion counter- 
parts. In fact their application is more cumbersome since there is no 
analogous recursion equation for the trigonometric functions. Therefore 

the first approach resulting in equations (2.3-1) and (2.3-2) is utilized 

in all but the most specialized cases since it was the first method to be 
successfully developed, and its application is more familiar than the sigma- 


omega form. 


EXAMPLE 1 (FEEDBACK COMPENSATION FOR LINEAR THIRD ORDER SYSTEM) 

At this point it will be advantageous to select a basic system for 
initial studies. The linear case is sufficiently unencumbered by mathe- 
matical manipulations to properly illustrate certain fundamental character- 
istics, and therefore serves as an appropriate introduction to the forth- 
coming discussions of section 4. Consider an uncompensated system with a 
forward transfer function of third order, whose gain is sufficiently large 
to make it unstable. Such a system is depicted by the unity feedback 
block diagram and characteristic equation shown in figure 2-1. The corres- 
ponding root locus diagram is included in figure 2-2 and it is seen by 
Routh's Criterion that the limit of stability is at K = 30,000. For 
larger gains than this critical value, it will be necessary to introduce 
either cascade or feedback compensation to regain system stability. Selec- 
tion of tachometer and acceleration feedback results in the compensated 
block diagram and characteristic equation shown in figure 2-3. Choosing 
a forward gain value of 100,000 results in a highly unstable uncompensated 
system. The stabilizing effect of either one of these forms of deriva- 
tive feedback without the other is shown by the cross hatched root reloca- 


tion zones of figure 2.4. The effects of simultaneous variation of both 
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forms cannot be as easily predicted on the root locus plane. At best, one 
can generate several root locus families, each as a function of only one 
variable parameter, and gain from these patterns an approximate idea of 

the required parameter values needed to place the dominant roots of the 
system in a prescribed position. The parameter plane technique becomes 

a very useful tool in this type of situation, however, since the pre- 
scribed root location automatically specifies the required value of zeta 
and also of W, . A cursory inspection of equations (2.3-1) through (2.3-4) 
will reveal that the following curves can be drawn in the parameter plane , 
which by convention has alpha on the abcissa and beta on the ordinate: 

1. Constant zeta curves, in which zeta is held fixed at a 
specified value, and w,, is allowed to vary through a pre- 
determined range. 

2. Constant undamped frequency or "natural frequency" ( W,, ) 
curves in which w, is held at a specified fixed value, 
and zeta is allowed to vary from 0 to l. 

3. Constant g= Wy curves, in which zeta and W,, are both 
varied such that their product is always constant at pre- 
selected curve values. 

4. Constant sigma ( o ) curves, which represent lines paral- 
lel to the imaginary axis in the complex frequency plane. 

5. Constant "damped frequency" ( Ww ) which represent lines 
parallel to the real axis in the complex frequency plane, 
and are constructed by holding Ww fixed, and varying sigma 
from zero through a pre-selected negative value range. 

6. Constant real root ( 7 +JO) curves, which represent 


specified fixed points on the negative real axis of the 
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complex frequency plane. (Thus it is seen that certain 
points on the s-plane are transformed into lines within 
the parameter plane. This factor plays a predominant role 
in singular line theory). 

Not all of these curves need be used by any means. The combination 
of any two of the first five curves positively identifies a complex root 
in the left half of the complex frequency plane, while the sixth type of 
curve locates the real roots on the negative sigma axis of the s-plane. 
Figure 2-5 is one type of parameter plane representation, in which Kg = 
alpha and K, = beta. It shows the constant zeta, constant W, and real 
root curves. One can see that for dominant roots at S$ = -16.5 + j35 the 
required values of alpha and beta will be at the intersection of the r = 
0.5 and Wy = 33 curves. In this example, they intersect at alpha = 
0.0005 and beta = 0.036. The third root is real and can be determined by 
the real root curve which coincides with these values of alpha and beta. 
It is seen that the real root is approximately at S = -80. These values 
are indicated by point A on figure 2-5. 

The curves for this example are deceivingly simple, and therefore 
easy to interpret. Additional care must be employed in higher order equa- 
tions, where all the root locations are not so easily defined. In order 
to ensure the desired dominant characteristics, one must certify the loca- 
tion of all existent roots. In exceedingly complex cases, this might neces- 
Ssitate a second, enlarged parameter plane plot centered at the point in 
question. This, however, is a greatly simplified task compared to single 


parameter techniques. 


2.4 Paramater Plane, Non-Linear Case. 


Equations (2.2-6) and (2.2-7) have been developed to handle two variables 
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when they are individual coefficients of a polynomial. Equations (2.3-1) 
and (2.3-2), or similar forms, apply to two variables embedded linearly 
in the polynomial equation. The following development is concerned with 
non-linear arrangements of two variables within the coefficients of a 
given polynomial. Consider first the case where the variables appear in 
product form. This arrangement is termed the alpha - beta product case, 
because the coefficients of the polynomial can be represented as 

ope ee CME SO ee +e. 


Substituting this into (2.1-5) and (2.1-6), one can write: 


2 K K z K K 
a 2 (No by Wy Uy, (8) ae 2!) cei a ee 
+ &B (1) “dy Li aitleeonpa Ge 


i= Cen Ww, A pe (¢) — O 
K=0 kK 
n 
Log, ene Oye tae ao es Cee 


nN) 


XB Z Ci dy Wy Uy CE) 


a Oe, ee Ce 


K=9 


Employing notational substitutions similar to those in the preceding sec- 


tion, we may write: 


Ex 2 CP Beemer, allo Cw 


0 K n 
2 : lo 
pee (1) bw, UD 
K=0 
n K e 
ai 2,09 Ce Wy, Uy €§) 
e / K 
Ca= ae ) Cy Ww), U , Ce 
Mac 


n 

ee ne oe ey 
n t 

D> = =. tall la) 
n Koo 

—— =, Cy WwW, nee C$) 

ee K 
E,= Ze Wr Uy (§) 


thereby reducing the preceding equations to the form 


AB, C8, Wy) + OC, (5, Hy) + BD, (5,4,) +E, (8%,) =O 


a B,(§,Wn) + BC2(S, Hn) +48 D2( 5) +E, (6,4) =O 


The functional subscripts have been included up to this point to emphasize 
the applicable dependent variables. It will now be necessary to elim- 
inate this descriptive subscript notation in order to achieve better 
comprehension of mathematical manipulations. Henceforth the dependent 
variables will be regarded as zeta and Wy, unless otherwise specified, 
and their inherent inclusion will be assumed in the remainder of this 


section. Accordingly, the above equations may be rewritten as 


CB me ECMO I 0 


a BEE Serre Be Dy SE, So 
The solution of these simultaneous equations for either alpha or beta is 


possible provided their Jacobian is not zero, where the Jacobian is de- 


noted: 


a (2 +6 D)(C, + ae = ac. See D,)\( Baker 8 
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Presuming this to be the case, the solution for beta is found as follows: 
fe, + @P,ILXB + 6C,+E,J =O 
[-B, -@BIi4 B,+ CC, +E, J=O0 


Le, D> en) 6° ae Ic, 2, -¢, By +E,D,-E, D, 1¢ ie FE B,-E,B =0 


which may be written in determinant notation as 


Bon P  (Rdeg +ALL IP F Ae, = O 
where 


A C, Cy Acg= C, C, 


D, Do | B, Re 
and so forth. Employing the quadratic formula, beta is found to be 


(Mea t Ag) ee | Ace Lies) * eee oe A ze) 


, Donbove 


a 
| 


(2.4.1) 


Proceeding in a similar fashion, it is possible to solve for alpha by 


the following process: 


ee, @ AD oh & But iGGpeot De), Cae FE, J =0 
[-¢,- oI Lek, +(C, tad e+ FE, T= O 


from which: 


(B,D, 7 B,D, oe” ie CB, C - B,C, Paeaale -E, D,) % es Bee = (0) 
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or in determinant notation: 
i _pn 


The desired solution for alpha then becomes 


Pe =. 2 ae 
ees |S (Asc +Aso) = (Ape * Aen)” -4( Asp Asc) =) 


2 Oep 





While the above solutions can be used to obtain alpha and beta, it is 
advantageous to simplify them for adaptation to computerized programming. 
Analyzing the terms under the radical signs, it is possible to show their 


equivalence as follows: 


\ (ge we = 4292) — (Dies + Aye)” utal timed 


Noting that Ac, == — A ec » and squaring both sides yields 


A pe +A,,)*- AA ep ee) =a Cia +A,,) — Gi Beg Acs) 


which, upon expansion, becomes 





(Age)? + Pee Ags) - Ca ~4 Dep Seo = 


(A ec a =e (Ape ied) a ee) ee 


Transposing and collecting terms 


4 Dep Oep Se ee ec + 4 Lep Age = O 


Dividing through by 4, and expanding determinants, 
Ges, Cy a By o Ce oS ayy, a (B,D, - B> D Ce) Co lege ) 


+ (C,D, -C, D,)(C£, 8B, -£,8:) = 0 


BC, OD, £, -8,C, Dye, 8, C7 ee, PBC), Eo 6, Gye eer 


+ B,C bE, —8,C, Dp, by + Bo G DE, — By Cope, 8 Gone) eee ee 
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Thus, we may define the terms under both radicals as Re ,» where 











\ 


Z j 
R= (Age + Ae) ~4HApp Oeil = ¥ (Aes tAgp)*-4Awde 


and the solutions for alpha and beta become 


a = pian he S 


B 





—_—_ oN BI aD) 
Y 


of the four possible combinations for alpha and beta, it is necessary to 
determine which of these satisfy the original simultaneous equations. 


The following results are obtained: 

WES & ¢ GeSoroDecuGe+ bale @ 

-D,£ B,x+ Ca BR + Dex +E, J = O 

(B,D, -5210,) 4 ec, Bh, -~0,D,) 6 + E,P,-E,D, = 0 


In determinant form this may be written: 


Substituting the above relations for alpha and beta 


(Asc tAes)* R (-Agc + Ago) +R , =O 
WE SE Rt Oe 
Z. "a 
Therefore, the radical coefficients for one term, either alpha or beta, 


must be the negative of the radical coefficient for the other term. This 


permits the following solutions for alpha and beta: 


» 
eee Ley ei Ne __ 
ee ee (2.4-3) 

- Bep 
(Aep — A ec) + (-1) BR ; 
ie. lake A Geia 1al ee Lr (2.4-4) 
P 2 Deo ' | 
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SINGULAR POINTS 

Inspection of (2.4-3 and (2.4-4) reveals that limitations are en- 
countered whenever A gp or pages equal zero. For these cases, it is 
necessary to resort to the original set of simultaneous equations for 
the complementary solution as follows: 

Should Agp be zero, the two solutions for beta are obtained from 
(2.4-4) and then substituted into the original simultaneous equations 


where the solution for alpha is obtained by the following relationship: 


iE 
oe nage? Ri a Sena ee f Cee 2 5 je (2.4-5) 
: By, D, + es Cy ot Bo 


When “\., equals zero, the analogous approach yields two solutions 
for alpha from (2.4-3) which are likewise substituted into the original 


equations, yielding the following values of beta: 





ee ay By + Ey _ %4BorEr le 
K we wrens 


= |,2 (2.4-6) 
Cyc, * L, CnC, +Ddz 


SINGULAR LINES 

The solution of a set of ™ simultaneous equations in Y} unknowns 
is not always obtainable. When a solution does indeed exist, it some- 
times possesses unusual characteristics in that it is not a unique solu- 
tion. These cases have been catagorized as singular lines, and their ap- 
plications to parameter plane theory are presented here for completeness. 
The following development is included for the purpose of showing the neces- 
sary and sufficient conditions for singular lines. 

Any set of simultaneous equations may be represented in the following 


matrix notation: 


[AJ LX] =LBIJ 
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where [A | is the coefficient matrix of order ™m x 

| X | is the variable matrix of order nx | 

| B | is the constant matrix of order mn xX | 

and [A B | is the augmented matrix. 

Any set of simultaneous equations is said to be consistent, (i.e., a 
solution is guaranteed to exist) when the rank of the coefficient matrix 
and the rank of the augmented matrix are equal. A unique solution may be 
found only when the rank of the matrix equals the number of unknowns, 7 
Otherwise the matrix is said to be singular. These three statements are 
very subtle, because their simplicity belies their importance. As an 
illustration of their application, consider the linear second order 


equations which were derived in the preceding section for coefficients of 


the form: 
Oe a en: ee (2.4-7) 








eum OnGs mlay y 
B, C “Ba 


N 


Thus for any characteristic equation of the form 

nN 

K 

ae Ae: a 

k =O 
where the coefficients are linear combinations of alpha and beta, it is 
possible to construct a planar representation as a function of the three 
terms, X, Y and ~ , with alpha and beta as the variable parameters. 
The inconsistency arises when the denominator determinant 2 is zero, 


and when at least one of the numerator determinants is not zero. Since 
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division of a finite number by Zero yields an infinite (or undefined) 
quotient, it is readily seen that the two simultaneous equations yield 
two parallel lines, which by definition are everywhere equidistant and 
therefore never intersect. Thus an inconsistent set of simultaneous 
equations has no solution. A unique solution is obtained whenever the ze 
term, or denominator coefficient is not zero. This is the most general 
type of solution, denoting the intersection of two lines at a finite 
point in the alpha - beta plane. Should one or both of the numerator 
deco nmiies X and Y be zero, it merely places the unique value on one 
or both axes respectively. The latter case of course is just a solution 
at the origin. The third classification for simultaneous equation solu- 
tions is the singular case which exists if and only if all three deter- 
minants X, \/, and ~ are identically zero. Such a case can exist 
when the number of unknowns in a set of equations exceeds its rank. Thus 
one or more of the equations is said to be linearly dependent. Geometri- 
cally this represents two collinear lines which by definition intersect 
everywhere, thereby giving an infinite number of solutions. Extensive 
value can be gained by the use of singular lines for design purposes, 
especially in self adaptive systems. It therefore becomes highly desir- 
able to be able to recognize the existence of such lines when depicting a 
system characteristic equation in the parameter plane. The theory of 
linear dependence Seascaatee that two linearly dependent equations differ 
only by a constant. Thus we may say that for (2.4-7) and (2.4-8) to be 


linearly dependent, the following conditions must hold: 


B, = KB, 
Co ey K F 
D 2S — kK Dy) 
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Because of the exceptionally involved calculations using the Chebyshev 
polynomial method or the binomial expansion method, the trigometric solu- 
tions (2.3-5) and (2.3-6) are amenable to singular line synthesis develop- 


ments. Employing those equations, (2.4-7) and (2.4-8) can be broken down 


into the following equivalences: 


n n 
roe t a eee ——- t ’ . 
ee Ke = Pell, «eciowie a. = oe = Se. sinté 
L=0 C=O : 
—_ n ( ‘ - nN l 
oC = kK a Se Cee ae z, Cp Sia 
—_ t= 
-_ yh ( n t 
KD, — KS d; Ay Coa o oo Db, = Z dj Wy Stnt@ 
(=0 


Using this criterion, it is possible to define the required coefficient 
values necessary to give singular line solutions, or to determine if 
any such combinations are indeed possible for a given polynomial order. 
Various extensions can be made to include all forms of non-linear 
polynomials. Such applications merely increase the number of trigono- 
metric equivalence equations, and involve conic sections rather than 
planar surfaces, provided the independent variables are limited to second 
order. For higher orders, no geometrical interpretation is possible; how- 


ever the fundamental theory remains unchanged. 


SINGULAR SOLUTIONS 


While most equations have explicit solutions, certain differential 
equations cannot be explicitly defined by specifying the values of the 
arbitrary constants. In Laplace or Fourier transform theory, this corres- 
ponds to non-parametric equations. The solutions to equations written in 


this form are called singular solutions [8]. Certain forthcoming studies 
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can be facilitated by expressing solutions for alpha and beta in non- 
parametric form. 


Up to this point, the simultaneous equations have been 
written as: 


xX, Ca, 6, en ee 


Xo (a, 7, 5 Wn) ae 


with their associated solutions for alpha and beta written as 
a= Ame 5. 225) 


B= BCS, wn) 
The following theory, 


therefore, is applicable to reducing the above equa- 
tions into the form: 


Vee Ct) meee 


with associated solutions 


GBs) 


and the alternate form 
ZC, #, Wr) ne 

with the solutions 

OSE Ce OD) 

6 = F(a, Wn) 


will now be shown that for all integers 


J and k 


z M3 


K K 
Cl) a, Wy U,., C$) = 0 (2.4-7) 

0 

whenever the relations (2.1-5) and (2.1-6) can be employed. 


That this is 
the case may be proved as follows: 


Given any polynomial of the form 
= Cl ae mes) mo 
mh possible to multiply both sides by ear c) f 
n K k 
SV l= Ve sGeicaalt rn eat tusiil §) Ue Eu 
k=o 


(2.4-8) 
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Given also that 


\3 


Ei) am, ae Ce =o 


a 


=0 


Multiply this by ys ( ¢) to obtain: 


~~ 


f 


Ee ’ nw Salen LU 
de a ir nN J 
K=C 


tlw ts) =O 


Subtracting (2.4-8) from this yields: 


K 


oe NN [ Uj (§) Ue) <9? = Us, (8) Un COMO 
K 


=0 


Recalling (2.2-3) we may rewrite the above in the following form: 

Z \ i (Ns | 

Se Cl) Sys L ee: Ca 70 (2.4-7) 
ie = 


With the aid of this relationship , it is possible to generate additional 
equations which will enable (2.1-5) and (2.1-6) to be written in non- 
parametric form. 

To derive non-parametric equations for constant zeta loci, it is neces- 
sary to eliminate 4), from the generalized simultaneous equations which 
then become 

Y Weckeg , ¢) = 0 
The necessary steps involve Sylvester's Dialytic Method of Elimination 

[9] for two polynomial equations of different order. An alternate approach, 
restricted to equations of equal order, can be used with Bezout's Method 

of Elimination but is not presented here [10]. Consider the following two 
polynomial equations of order ™ and % where ™ may not necessarily 
be equal to ‘vi , and where the coefficients of omega are functions of 


alpha, beta, and zeta. 
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a 2 Ta =| ce 
7 ( WwW, ) = oe ta ae Wan + ro. Wn + 2 8s 8 ee WwW, ~ i w, 
yn 
a oo 
= z Fy (% 8,4) Wy,” =O 
== ve 4 N- 
Yo (Wy) = oe og | Wy i 92) oe a eee +9, wy," 


2 ik 
= 2 9x 0%, 8) §) Wy = O 

K>o0 
It will be necessary to form a square coefficient matrix out of the 
combined equations to enable the determinant to be taken. This deter- 
minant is then set equal to zero in order to obtain the non-parametric 
solution for all cases but the trival one which occurs when W, _ simply 
reduces to zero. The square matrix is of order mxm_ and is formed 


by writing vi (W,) and Vo (W,) as follows: 


=| 
j n~| m kK . 
vs (w,) = S Wy LY, (wy)] A S F (8, §), W,, Ww,” = OM 428) 
J=0 g20 =O 
m=} ij ames ae n , 
V4, (Wy) =H win | Ye (wp)| = = 9. (% @, 9) wl Ww, =(0 (2.4-9) 
J=0 i at 


2 
This corresponds to multiplying each term in Y , (4,) by i Wry) Wy ee 


yin Nae Bay 


W, ) Yn » Hn respectively. 
In a similar fashion each term in ,(W,,) is multiplied by !,“n) 
2 -! 4 ; , Sree 
Wry sree Ww respectively. The matrix equivalent is shown for amplifi- 


cation on the next page. 

Equations Y (On and Vo (Wn) were not constrained to be of equal 
order, so that the preceding steps might show the method to be appli- 
cable for unequal cases. Practical use of this procedure for the elimina- 


tion of Wy, in context with the parameter plane, however, results in two 
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TT 
0 O o i i f, r ; Sy ~2 fin: } a | W . 
CeO Geet ie tar ims! Ime eC | ww! 
O ft F, Sy + fe Giete oO Oo | | Tg 
re 
% 7, 3 : : Fin ‘wae Thee Oc O | | Ww? 
6 0 60 0 0 6 Bs, G 5s | | ws 
na) = ‘ | wag 
y 6 0 0 VOU? i “Tike ° | = 
2.4-10 
So ome ft ho — 
| 
2 a a a ee aC aa 
| 
PO ma a” ER Ce Oe One 
g. 8, J3 : 4», G G 6 OC 0 | 
Jo ap Jy ‘ ae c C 0 eo. 06°0 0 | oo 
a 


— 


equations of equal order when the following substitutions are made far f,, 


and qk ‘ 
me Gl) a, u, (5) 
3 = (No, Pons (§) 


The non-parametric equation for constant W, is derived by tactics 


analogous to those used to generate the constant zeta equation. 


Again, 


a set of equations is obtained for the purpose of taking the determinant 


of the coefficient matrix which, by definition, must be a square matrix 


if the determinant is to exist. Since the variable matrix U .(§) is 
only of the first power, it is possible to generate Nn-~! additional 
equations by use of (2.4-7) which was previously derived. Letting j 


range from zero to Y-/ in (2.4-7), the following set of simultaneous 


equations is obtained: 


a 
= 
x 
9 
a 


Wy, ye) ao 


9 


Ei" &y Wy” Ux) C$) = © 


Ss 


Se 7S 


(-1)" a. Oe U ng (§) a 


Ss 


EC) a, Wy" Un l =o 


= 
" 
© 


a), Oe ae oO 


x 
o 


(2.4-11) 


2 (-1)" ay Wy " U ka(n=3) Lo 
z é1)" i Wr Ue -(n-2) eo 
Z é))" i il Uk -(n-1) (S]= 0 


Writing these equations in expanded form: 
CC ine Qo LO U, C§) at. Ci) a, WU, UG (¢) + Cine, We La) Flee Li) "an Wy Up (6)=0 
Gail). emiblamml te Ge luau itl O i tenet) aoe), (emGg) cciemmes . ba, W, Un =O 


(-1)°A, Wy’ Ung (G) + EN) 4%, Wy Uy (8) + C1)", Wy Uo (€) + 22 CDW," U9 (9 =0 


(envi) Umea MN), ng i Dag Re 


. +(-) a, wow," U, Of) =O 


44 


These equations may be simplified by substituting the following 


relations; 
Lo (¢) = 
U_, () =-U, 


Thus we may write: 
n n \ 

O + 1) Ct, Ce) ee Ce eet re a 
mC iy Og ve ul Cs) ae O anagne 7) Gaon ere 
—(-1)" a, Wy? U2 (4) soba. jmwelyemel Gale joome ~ommnde (-1)" ay Wy Uqr G)= 0 

e | n n . 
nif =!) Qo tig Cail diel L) me (-1) Qa, Wy Maceo) I aa (-/) an Wy, U, (§) =O 


Which is seen to be a set of n-1l equations with n-1l unknowns. In 


matrix notation these become: 


. J Zz 7. n ” 
(EG Uy ae (-l) AgWy tit’ Cl) ay, Yp U, 
I 2 ci 3 
C1)" a, Wy, — f-!) Qo Wr, i Cie a3 Wy fawense U> 
| 3 ! i a ie? ee ae 
Ee a, we —(-1) @ Wy 4 (-1) Ay p= (I) Igy fever U, 
= 

a 4 ze Zz 5 5 ‘ ' 
C1) 24,W, -f-) a w + (-1) 4-Ww, ~(-1)O,Uq +00 ble 
i pre yecey 2. Wn on Pia) 


n-2 n-2 


(1) an Ya) Oy os 
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ee EE 
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- -----+ + 


?-v 
= 
Ss 
> 
Ss 
e 
s~ 
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th order linear- 


The above developments have shown how a set of two 1) 
ly independent equations with four variables can be reduced to one equa- 
tion of order n-l with only three variables. This is accomplished by 
solving the coefficient matrices of (2.4-10) and (2.4-12) which elimin- 
ate zeta and omega respectively. The application of these methods will be 


indicated in the following article where the mathematical concepts of the 


frequency response curves are presented. 


2.5 Frequency Response Applications of Parameter Plane Theory 

Frequency response curves can be portrayed in many ways. Some of the 
more useful ones will be considered in sections 3 and 4. The various 
representations of frequency response curves can be derived from the 
following fundamental theory relating magnitude to frequency. Consider 
again the transfer function of a mathematical model which represents the 
system to be studied. This can be reduced to a ratio of two polynomials 


by partitioning, or other suitable means to arrive at: 





Nn . 
T 6s) = eS NO) oy 
Ea +b, S + bgs* + bg site by J Dats) 
eis 
Jeo 


To find the absolute magnitude of this transfer function, the following 
is noted: 

ZL . re 
= T iwi] = [ Fi) | TT -iw)| 
(Although the notation standards of section 2 clearly state the distinc- 
tion between W and Wy, , it is important to note that natural 
angular frequency is no longer being used.) The following steps will 


illustrate this transition: 
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i) Fg ee 


\2 


a a E TCL + dies GUC) + a {jw)? 4 a ky (jiv)” 





Rewriting this to separate the real and imaginary parts, we have: 


‘ 2 ‘ 1 ‘ 3 
Qe ta, (see) a (fe + Qa,jw +a, (ju) po 





ow): = hint...” null. an cn 
be +b, CO + by, (Ww) eee? Sie als CI. 
[ a. -a, w* +a, w"- pie] eo abe ao” teal 
TQiw) = =— a oe 
|b. -b,w byw J+ jf bw bw? +e 


In abbreviated notation, this becomes: 
_Neven (:¥) + Nodd (iW) 


Dive we). midi Lai) 
Proceding in a similar manner we find the following: 





or 


Qa, -a Pee ue —"- : a Q,W-A,W t+ aw, 
T (-jw) = [2s - a. eoiis Z anes elin 








or £ 4 
b= BW + byte +} Lew byw? + be we 


Employing abbreviated notation again, we may write: 


Neven (iw) —Nedd (iw) 
T (-jw) = SS 
Deven (iw) — boas (7) 
Thus the square of the absolute magnitude is found to be: 


[Neven (ju) T Nodetio) |i Nevertr one ccd iam 
[Dever (iwi + Dodd Cw] f Deven (iw) — Dodd Ciw)] 


2 ee (jw) ]* om [ Nodd ‘alee 


ee ae (2.5-2) 
[ Deven Gw)}” =» | Dedd(iw)]* 
& Beo- 4.u* +a, ieee | "Taw a,W? +a. wre)” 
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This can be written in the following summation notation; 


Zn K k Zn-| KS K 
|) Teme + = (-I)~ ¢,w 
= K=/ 
Me — Weve: K odd 
—_ 2m £ > 2m = | a | 
& g ea L 
ae 2 
ey Cl) * daw 1 . ener io 
Leven ad 
n 2k 
. Z (-1) Con W 





Let 

2 
X= Ww 
eK = Cok 
Ts = dog 


Substituting these new variables above, a more concise notation may 


be formalized for easier illustration of pertinent relationships. 


K K 
a 4 
z C}) Cb. aX 


M” — -kee 


2 L 
= (—1) Ff, x 


(2.5-3) 


In general, the coefficients of (2.5-1) are non-linear combina- 
tions of at least two variable parameters, designated alpha and beta. 
The linear combinations were discussed in article 2.3 and the associated 
form of any polynomial coefficient was found to be 


— 


K by o ro Goer —+- any 


Q 


For the non-linear case, the simplest combination was seen to be the 
alpha-beta product combination. This subject was treated in article 2.4 


where the resultant coefficients of a given polynomial were broken down in 


the following manner: 
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On= byt + Ce B + Ge AR + eK 


Theoretically, there is no limit to the degree of difficulty which exists 
in the non-linear case. A few examples of the complexity involved in 


non-linear parameter combinations are presented here for completeness. 


aa ial a a i a a a 
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Clearly, it can be seen that such combinations are exceedingly unwieldy. 
Attempts to handle all but the simplest arrangements are extremely labori- 
ous. Computer programming is still in the embrionic stage for parameter 
plane applications. Attempts to utilize this tool have been extended 

only to the alpha-beta non-linear case. This restriction is acute when 
considering the frequency response curves. A cursory review of the 
mathematical developments at the beginning of this article will show that 
the polynomial coefficients of (2.5-3) are obtained by squaring the poly- 


nomial coefficients of (2.5-1). Thus, even the simple linear case of 
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article 2.3 results in non-linear combinations of second order for fre- 
quency response studies. To handle this problem, it becomes necessary 

to depart from the procedures which were employed in previous develop- 
ments. By assigning a value to either alpha or beta, the task is reduced 
to merely a quadratic formula solution for the remaining unknown variable. 
This can be illustrated in the following manner: Given a non-linear co- 
efficient combination of the alpha - beta product form: 

Oy bet. eG ee ecyeeceeeet (e) 

Assign a value to alpha, and denote the value as A: 

ee b, H tC Geo eh ee 1 ey. 


This may be rearranged for clarity, 


Qa [ bpA + ey ae fce,+d,A] @ (2.5-4) 


where the first term on the right hand side denotes a constant value, and 
the second bracketed term represents a coeficient of the independent vari- 
able beta. In a similar manner, a fixed value may be assigned to beta. 
For illustration purposes, let it be denoted as B. This gives: 

Co ec ene Cy See ey 

which is likewise rewritten as: | 

Q, = Le B + e,] + [bye + dB] % (2.5-5) 
Again, the first bracketed term represents a constant while the second 
denotes the coefficient of the independent variable alpha. The computer 
program for frequency response oe program PARAMS, is based upon this 
solution procedure. Consider first the case where magnitude and omega are 
specified. By assigning incremental values to either alpha or beta, the 


corresponding solution may be obtained for the other. This is the method 


employed in the graphical representations which have both alpha and beta 
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on the axes. For the graphs where magnitude is an axis variable, it is 
incremented while a corresponding value of beta is obtained as a function 
of alpha and omega. The same philosophy is utilized when omega is an 

axis variable. For the case where both magnitude and omega are on the 
axes, the values of alpha and beta are pre-selected and the familiar Bode 
plot results. All of these graphical arrangements are treated with con- 
siderable detail in section 3. Their introduction at this point was 
simply to point out that the solutions are based on (2.5-3) and either 
(2.5-4) or (2.5-5). Since (2.5-3) is the fundamental gain versus fre- 
quency equation, it can be anticipated that all frequency response para- 
meter plots will contain the same basic information as a series of con- 
ventional Bode frequency response curves. When two variable parameters 
alpha and beta are utilized together with magnitude and omega, the problem 
becomes one of four dimensional space. When still another polynomial co- 
efficient variable parameter is considered, an additional dimension is 
added. In spite of this, it is still possible to derive one set of para- 
meter curves from another, since each is just a different projection from 
the same multi-dimensional model. The value of each separate representa- 
tion lies in the ability to depict certain important characteristics which 
would not be obvious otherwise. Additionally, certain configurations 
facilitate the study of a dominant parameter, or the effect of certain 
parameters on different performance requirements. Certain mathematical 
considerations must be pointed out at this time concerning the inter- 
dependence of these graphical plots. Clearly, the construction of PARAM-7, 
the Bode plot, offers no problem. In this case, the desired values of 
alpha, beta, and omega are inserted into (2.5-3) and the corresponding 
magnitude squared value is then computed. For the remaining programs 
(PARAM-1) through (PARAM-6), a rearrangement of (2.5-3) is necessary since 


aL 


magnitude is no longer the unknown variable. As an illustration, con- 
sider (2.5-4) where alpha has been specified. Beta then becomes the un- 
known value to be solved. For added clarity it is better to revert back 
to (2.5-2) in order to point out why a second order system is always en- 
countered in the linear,case, as well as the non-linear alpha - beta pro- 
duct case. Since beta is the unknown variable, pre-selected values of 
magnitude, omega, and alpha must be presumed. Performing these substitu- 
tions, (2.5-2) may be written as 
yr Ltt nO)” +0 % +5. 6)" 

ellie, lige [vy +m Pp)” 


where, using the notation of (2.5-4) 


(ee | bow or Gy ee baw es we + [ b, A rege = cs 
Gece POT = 1, + ao Aw fc, +d, A} wi -- 
Tere, ee, Ate, | Pe Lb, P ve 


| 
| 
| 
| 
2 4 
aes Pare ee a Atw + Pep eal aoe 


Therefore, one can easily see that by squaring the right hand bracketed 


values and combining like terms, 


Mra SNBO +SNBI 6 + SNB2 p* 

a De aS 8 15.049 0) 8? 

where the terms SNBO, SNBI, SNB2, SDBO, SDB1, and SDB2 are so chosen be- 
cause they are the variables used in the computer program to be dis- 
cussed later. Since the numerical value of mM? has been specified, this 


equation can be cross-multiplied. Transposing the right hand side of the 


resultant equation, one can write: 
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SN8O ~M* SoBO +[ SNBI-M* SpBIT@ +[5NB2-M*SDe2] 6'=0 
This is clearly a quadratic in beta, and employment of the quadratic equa- 
tion yields the desired results. Imaginary numerical values have no 
physical significance for graphical plotting and are therefore discarded. 
However both the plus and minus signs in front of the square root term 

of the quadratic formula must be considered. Thus it is possible to 

have two values of beta for a specified combination of magnitude, omega, 
and alpha. This is an important consideration in graphical plotting 
techniques which can lead to spurious lines if not properly accounted 
for. Effectively the problem requires calculating all the prescribed 
graph points using first one sign of the square root term, and plotting 
them. Then the computational process is repeated with the alternate 

sign of the square root term. This illustration has shown how a solu- 
tion would be obtained for beta when alpha, magnitude, and omega are 
known. Precisely the same procedure would be followed to obtain a solu- 
tion for alpha when the values of beta, magnitude, and omega are given; 


only (2.5-5) would be used in lieu of (2.5-4) 


RESONANT PEAK MAGNITUDE 

Performance specifications of certain systems include stringent re- 
quirements on the presence of a stipulated resonant peak magnitude. Con- 
versely, other design studies are concerned with the guaranteed absence 
of a resonant peak magnitude above a pre-determined level. Either of 
these requirements results in the study of frequency response curves with 
particular emphasis on the peak magnitude and the frequency at which it 
occurs. To facilitate such studies, it would be advantageous to construct 


resonant peak magnitude curves for the system as a function of the 


53 


polynomial variable coefficients, alpha and beta. This involves the 
utilization of the non-parametric theory presented at the end of article 
2.4. 

Normally, one thinks of only one maximum point in the entire fre- 
quency spectrum because of dominant root characteristics of low order 
polynomials. Certain high order systems, however, can be represented by 
characteristic equations which have more than one pair of dominant com- 
plex roots. Such arrangements will yield more than one maximum point on 
the frequency response curve. Only one of these points can be consider- 
ed as the global maximum, while all other peaks must be relegated to the 
role of local maxima. To avoid ambiguities, the notation M p (wv) will be 
restricted to the global maximum point only. The corresponding value 
of omega will be denoted by W > . Since it is desirable to refer to 
local maximum points as well, the term M » (w) shall be construed to be 
applicable to all resonance points, both local and global. The respec- 
tive angular frequencies at which these occur will be labelled as UW 
Thus it is possible to have any number of My (w) and WU), values for a 
given transfer function. Only one of these can be the M p@) and Wp ; 
however. 

Recalling (2.5-3) where aN and on are both functions of the vari- 
able parameters alpha and beta, it is possible to state this equation in 
the following functional notation: 
fF, (a, 8, Ww) 


M*(w) = 
§, (a, 6, w) 


(2.5-6) 
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Assuming at least one pair of dominant complex conjugate roots, it is 


possible to find the following relation from (2.5-6): 





M a = - faa (o, 6, WW) aa vespu ay ioe (4, F, wi) ae’ 
Fo (ah, 6, w) Fy (a, By wy) 
W= Wy 
a (2.5-8) 
70, 6, 2) f¢ Col, fos We) 
W=Wp 


Fundamental theory of equations dictates the following correlation be- 
tween (2.5-6), (2.5-7) and (2.5-8): 

1. A relation in the form of (2.5-7) exists whenever the 
slope of (2.5-6) is equal to zero and its radius of curvature is nega- 


tive. Expressing this statement in mathematical terms: 


a[m* cw) | 


_ 4 (2.5-9) 
gw =e & 
2 
\ lua 
aL mM a (2.5-10) 
o 0 


2. Only that value of Wy, in (2.5-7) which globally maxi- 


mizes (2.5-6) satisfies (2.5-8). Thus it is possible to state: 


( uw ) oa MAX imum ; M> (w i (2.5-11) 


) 

ae 
The preceding theory is adequate in that it gives the necessary and suf- 
ficient conditions for the existence of M p(w). Satisfying the above 
criteria is fairly straight forward whenever the values of alpha and beta 
have been selected or prescribed. The dilemma arises when seeking the 


2. 
algebraic relationship for Mp (w) as a function of alpha or beta. 
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Although three equations are available, the inequality in the second de- 
rivative criterion of (2.5-10) imposes a severe limitation on the solu- 
tion technique. The most straight forward approach would be to utilize 
(2.5-7). From these it is possible to obtain (2.5-11). Such a method 

of attack is only useful when both alpha and beta are given. This would 
require a two dimensional scan for all possible values of alpha and beta. 
Various computerized programs exist in optimal control theory for the 
determination of (2.5-7). These, however, are based on apriori knowledge 
or random search techniques using Lagrange multipliers. Not only are 

these methods time consuming, they are not conducive to highly sensitive 
frequency response curves where resonant peaks are in close proximity to 
one another. An alternative solution employs the non-parametric theory 
which results in (2.4-10). This procedure is applicable only when the 

two equations involved can be set equal to zero. Therefore it is neces- 
Sary to investigate the locus of inflection points of the magnitude versus 
frequency curve. From this information, it is possible to determine how 
many local maximum points exist. It can be shown that for 7” inflection 
points, there can be no more than % #1.maxima. In practice, however, only 
{rom will be found for most systems. This procedure has the important 
restriction of eliminating omega. Thus it is impossible to employ (2.5-6) 
in order to find the value of ™ *¢w) at each inflection point. At best, 
this result makes it possible to construct the locus of inflection points 
in the parameter plane when alpha and beta are the axis variables. It 

then becomes necessary to employ some other means of determining which side 
of this locus represents the minima, or maxima. Full utilization of the 
non-parametric theory can produce this desired result. Inspection of the 


developments leading to the coefficient matrix of (2.4-12) will show that 


56 


the combined application of (2.4-10) and (2.4-12) provides the necessary 


information for (2.5-6) as follows: 


l. The results of (2.4-10) yield a solution of the form: 


= eee (2.5-12) 


2. The results of (2.4-12) provide the following alternate 


solutions: 


y = ae ( we i) (2.5-13) 


/ 


or 


| 


£ , ( é, Ww ) (2.5-14) 


depending on which coefficient variable is eliminated. From these re- 
lations, it is possible to employ (2.5-6) directly by use of (2.5-12) 
and (2.5-13), or (2.5-12) and (2.5-14). In this manner, it would be 
possible to employ the optimal control theory previously mentioned. 
However the more direct approach would be the utilization of the gradi- 
ent technique which is presently employed in many optimal control pro- 
blems. Alternately, either of these pairs may be solved simultaneously 
to provide a locus of inflection points as a function of the two re- 
maining variables, alpha and omega, or beta and omega. This latter ap- 
proach would be most desirable for frequency response parameter plots 
employing these variables as the graphical axes. 

The preceding theory has shown several ways of handling the peak 
resonant magnitude problem. Inherent in each of these, however, is a 
vast array of mathematical manipulations which is limited by solution 


techniques presently available. The following development will indicate 
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this more clearly. Consider again the fundamental relationship of (2.5-3) 
upon which these derivations are based. 


The following result is obtained 
for (2.5-9): 
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Let Q;= 


numerator coefficient of (2.5-15) with the sign embedded, and 


let hg = denominator coefficient. Thus (2.5-15) can be rewritten as: 
yn ey} 


dM? (x) agen at aieals 
o) x 





m A 24 
é (-1) me x 
The equivalent form of (2.5-10) then becomes: 
m a 2L ada . a 
a2: - 
Be Mae (x) | 200 h ya ae Tat Kom) 


ca aciital slieoedaaiialen 


(2.5-16) 
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Thus from (2.5-16) it is noted that the order of the denominator poly- 


nomial is raised to a power of four times the order of the denominator in 


the original transfer function. In like manner, it is evident that the 
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numerator polynomial of (2.5-16) is also of very high order. Therefore, 
all but the simplest of transfer functions would require computer assist- 
ance to solve the determinants of (2.5-15) and (2.5-16). Although many 
programs exist for explicit determinant solutions, these cannot be used 
when one of the terms is an unknown variable. Until a generalized com- 
puter program is developed to handle these cases, the construction of peak 
resonant curves can only ba achieved by a laborious calculation of second- 
ary curves. Such procedures do exist, whereby tangent lines are drawn 
from a series of curves. To be highly accurate, however, an excessive 


number of these are required for sensitive systems. 


PHASE CALCULATIONS 

The preceding discussions in this article have been limited to magni- 
tude studies in the frequency domain. An equally large field of interest 
is that of phase angle characteristics. The utility of the magnitude 
versus frequency plot lies not only in analyzing or predicting the fre- 
quency response of a system, but also in the ability to check its stabil- 
ity as well. The correlation between magnitude and phase curves provides 
this verification by means of phase margin and gain margin calculations. 
It therefore becomes desirable to derive the phase calculations as a func- 
tion of the three variables: alpha, beta and omega. 

Initial developments in this article showed how (2.5-1) could be alter- 


ed to produce 


\j Eien Gi a) aa odd (jw) 


[ac] ee 
bi even (J4) + Dodd (jw) 
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Considering the numerator and denominator separately, either can be re- 
garded as the real and imaginary components of a vector with a given magni- 
tude and angle. Magnitude calculations have already been developed and 
are of no interest here... Phase angle calculations, however, are known 


to involve the arc-tangent of the vector components as follows: 





Phase Pngle — QA = fan (2.5-17) 





Re { TU) $ 


Two solution procedures are therefore possible. One is to solve independ- 
ently for the angle of the numerator vector and subtract the denominator 


vector angle from it. Thus: 


= ie lll at rite Al 


0 = tan |_Netd Cw) “Y) Dodd Civ) 
wan! 2 cele) ua) Deven (jw) (2.5-18) 


The alternate approach involves only one trigonometric conversion. Clear- 


ing the imaginary term from the denominator of | (jw), 


Neven(iw) + Nodd (3w)| | Deven (iw) —Dedd Civ) 
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Therefore, substituting this into (2.5-17), 






=I ae (jw) Deven( jv) — Neven (jw) Dodd (iw) 
Neven(jw) Deven(iw) — Nodd (iw) Dedd(iv) 





= fan (2.5-19) 


Either (2.5-18) or (2.5-19) may be used for the prescribed phase calcu- 
lation. For the Bode plot (PARAM-7), both have the same degree of diffi- 
culty in computer applications, although (2.5-18) is preferred on long 
hand calculations. Extensions to phase curves in the remaining para- 
meter plane frequency response graphs (PRAAM-1 through PARAM-6) have not 
been attempted. These would be possible without the use of the quadratic 
formula if (2.5-18) were used. Selection of (2.5-19) requires the usage 
of the quadratic formula with both its square root coefficients as il- 


lustrated in the magnitude portion of this article. 
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3. Computerized Graphical Construction Techniques For Frequency 
Response Studies. 


Graphical analysis has always held a virtuous A dindist on in all 
facets of industry and business. Since the Second World War, however, 
the advances of technology promoted the engineering profession beyond 
the existent graphical methods which were derived for the single purpose 
of varying only one parameter at atime. In spite of this, the ancient 
adage "'one picture is worth a thousand words" hag withstood the test of 
time, as witnessed by many practicing engineers today. Many pictorial 
representations of single variable systems are familiar; however with the 
advent of multiple parameter variations, graphical representation in two 
dimensions becomes more complex. Considerable experience and moderate 
imagination are necessary to visualize three dimensional configurations 
when presented in two dimensional planes, in order to adequately make an 
intelligent intrepretation for meaningful results. 

The recent developments employing algebraic methods for variation 
of two parameters have resulted in the successful use of digital com- 
puter facilities to perform the required manipulations for graphical 
representation of mathematical equations. The preceding section was 
devoted to a fundamental development of the underlying mathematical theory 
employed in two parameter design studies. This section will include a 
broad presentation of certain graphical methods which are most likely to 
be employed in two parameter electrical control problems, with primary 
application to frequency response studies. This section is also con- 
cerned with a detailed inspection and description ef ene computer pro- 
grams utilized in formulating those graphs. The discussion is not direct- 
ed toward programming, but rather an effort is made to clarify certain 


computer oriented problems so that engineers who desire to utilize any 
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program contained herein may do so with minimum difficulty and yet gain 
maximum benefit from its comprehensive construction. Illustrative ex- 
amples of each of these programs appear in sections 4 and 5 where several 


case studies are investigated in depth. 


3.1 Selection of Graphical Representations. 

For any given set of six variables, the study of their respective 
inter-relationship by means of two dimensional projections results in a 
myriad of graphical portrayals. The magnitude of the problem can best 
be seen in the following development: Assume that an arbitrary number 
of variable parameters are to be considered. If 7 is the number of 
design variables, it is possible to show that the number of possible 
graphical representations is M-factorial ( 7 ). All graphs reveal 
essentially the same information content when one axis variable is switch 
ed from the ordinate to the abcissa, and vica versa. Therefore, redun- 
dancy is eliminated by allowing the abcissa and ordinate variables to be 
interchangeable. This reduces the number of combinations by a factor of 
two, thereby making the revised number of graphical representations equal 
pow: ¢ nal ). Although a considerable reduction in the size of the ar- 
ray has been effected by interchanging the axis variables, it is highly 
desirable to reduce this size further, since a six variable array means 
360 graphs, a five variable array requires 60 representations, and so 
forth. To further reduce this compilation, it is often possible to inter- 
change any two parameters alpha ( “ ) and beta ( © ). Such would be 
the case in the parameter plane technique where they are, by definition, 
arbitrary variables of the characteristic equation. Even this reduction 
leaves a formidable task of evaluating over 180 graphs for a system in- 


corporating six variables, and any system with five variables requires 


68 


an array of 33 different graphical representations. For this reason 
alone, the present state of the art is limited to at most four variables 
for practical design techniques employing graphical displays. Any system 
with four variables requires basically seven different graphical presenta- 
tions for comprehensive evaluation. Seldom are all of these required for 
a satisfactory design, although barring the elements of time and expense 
they would most likely be necessary for optimal results. Effectively the 
problem at hand involves the inspection of four independent variables. 
This can be regarded as a four dimensional space which is not physically 
realizable. By assigning a constant value to one of these four variables, 
the problem is reduced to a three dimensional sub-space representation 
which does have a physical shape. It's meaning, however is restricted to 
that particular value of the fourth variable. To present this three di- 
mensional sub-space on a two dimensional plane requires that a large 
number of curves be projected onto this plane, where each curve represents 
a different value of the third variable. The following discussion is 
based on this concept, and the resultant graphs should be interpreted ac- 
cordingly. 

Most specifications for the dynamic performance of linear systems are 
written in terms of the time domain and usually include such terms as peak 
overshoot, transient duration, and steady state error. It therefore be- 
comes necessary for the control engineer to bridge the time domain with 
the more utilitarian frequency domain. Correlation is accomplished by 
means of the Fourier integral when rigorous methods are called for. This 
treatment is avoided in practice however, because of the inherent arduous 
task involved in all but the simplest of integrations. It can be shown 
that by suitable approximation, many linear problems fall into the cate- 


gory of second order systems. For engineers not experienced in partitioning 
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methods and root locus techniques, this statement would seem to be easily 
contested; however, such an analogy actually exists in most cases. For 
systems which cannot be resolved into an equivalent second order approxi- 
mation, the absence of dominant complex roots is evidenced by the shape of 
the closed loop frequency response curve. Since considerable information 
exists in present day literature on the panerati zed second order system, 
much time and labor is saved by manipulating a given problem into this 
form whenever possible, and extracting the equivalent damping and band- 
width specifications from pre-computed curves and tables. 

A cursory review of the terms used in control engineering applica- 
tions will reveal that some of the more common design criteria are steady 
state magnitude ( M_ ), phase angle designated theta ( ©  ), natural 
angular frequency designated omega ( Wy, ), and damping ratio designa- 
ted zeta ( Some problems are more amenable to the rectangular 
coondl nates of Che wosnnter frequency plane. These are designated as the 
damping constant sigma ( co ) and the damped angular frequency omega 
( oo ). It is again emphasized that two forms of angular frequency are 
designated by the word omega. Much of the existing literature uses this 
term very loosely, and extreme care must be exercised not to mistake un- 
damped, or natural angular frequency ( WW, ) for damped angular frequency 
( W ). The relation between these two variables is found in article 2.1, 
and it should be noted that these terms are synonomous only when the damp- 
ing ratio ( ¢ ) equals zero. 

Section 2 developed the notational concept of alpha ( o& ) and beta 
( @ ) as representing two arbitrary variables embedded within the co- 
efficients of the characteristic equation. The preceding paragraph has 
shown that eight variables are frequently encountered in the resultant 


electrical engineering control problems. As previously stressed, the 
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simultaneous usage of all of these is not feasible. Fortunately, however, 
one may break the control problem down into two basic categories: stabil- 
ity and frequency response. These two facets are interrelated, however, 
and for that reason some attention must be given to system stability 
when performing frequency response studies. 
The parameter plane stability problem may effectively be considered 

as the study of the relation between the characteristic equation poly- 
nomial and the stability performance variables { and Wy, , or @ and 

J . When the two parameters alpha and beta are embedded within the 
characteristic equation, the scope of the problem increases to a study of 
at least four variables. Thus one might be inclined to consider seven 
possible graphical arrangements. Unlike the frequency response studies, 
however, the stability problem has been attacked from a slightly differ- 
ent viewpoint. Only alpha and beta have been plotted along the axes, and 
no attempt has been made to develop the necessary equations and programs 
to portray a different set of graph coordinates. Thus only one basic 
graphical presentation is used for parameter plane stability studies. The 
example in section 2 illustrated this method of stability study by employ- 
ing the following alpha - beta parameter plane curves: constant zeta ( ¢ ) 
curves; constant omega ( “y ) curves; constant zeta - omega ( [-W,) 
curves; constant sigma ( OG ) curves; constant omega ( WU ) curves; and 
real reet.( a + j ©O ) curves. Although the constant zeta - omega (§-Wp) 
curves are identical with the constant sigma ( @ ) curves, the alge- 
braic methods required for each are different. The names assigned to each 
have therefore been kept separate to avoid confusion in the respective 
graphical construction techniques. Additionally, the real root (o +) O ) 
curves represent constant sigma points in the complex frequency plane which 


are transformed into lines on the parameter plane. As such, they are not 
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to be confused with the constant sigma ( @ ) or constant zeta - omega 

( : - Wy, ) curves. Although partial programs exist for these curves, a 
program entitled PARAM-D has been written for the primary purpose of 
consolidating these programs into one package for added flexability and 
versatility. This permits any combination of the above curves to be pro- 
cessed on one graph. 

Complex engineering problems often involve a stability study of systems 
with more than two arbitrary characteristic equation coefficients. Pro- 
gram PARAM-D has been extended to include a third variable which is here- 
in defined as gamma ( Y ). This scheme will allow stability analysis 
studies to be visualized in three dimensions, rather than two. When con- 
fronted with this type of problem, it then becomes essentially a four 
variable problem. Therefore it is necessary to hold one of these four 
independent variables constant, and generate a series of curves using 
the remaining three. This is accomplished by assigning a constant value 
to one of the stability performance variables: ¢ Wn) f-w,. 7, Ww, or o+j0. 
The remaining three variables are those in the characteristic equation co- 
efficients, namely alpha, beta and gamma. With these, the program at- 
temps to develop a three dimensional surface for the fixed value of the 
specified stability performance variable. The x, y and z axis variables 
are alpha, beta and gamma respectively. Graphical presentations are re- 
stricted to two dimensions, however. The dilemma therefore is that this 
three dimensional surface must be projected onto a two dimensional plane, 
which is herein defined by the x and y axes. The only way that this pro- 
jection can be accomplished is to construct a large number of third vari- 
able (gamma) curves on the graph. Each curve then represents a cross 
section of the three dimensional surface. The physical shape of the sur- 


face can then be visualized providing a sufficiently large number of 
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these curves are thus constructed. 

The frequency response study has been viewed in a manner which is 
slightly different than the above stability program approach. It was 
initially decided to make a thorough analysis of the four variables: X , 

Br > M , and W ._ As such, alpha and beta can be regarded as two 
variables embedded in the characteristic equation coefficients, while 
magnitude and omega are the two frequency response performance criteria 
of interest. Treating these as four independent variables results in 
seven basic graphical combinations. These have been titled PARAM-1 
through PARAM-7. Additionally, phase is investigated in PARAM-7, where 
its correlation with magnitude is most easily made. All possible graphi- 
cal presentations of interest have been utilized, and a summary of these 
is presented in figure 3-1. They have been implemented in one consoli- 
dated computer program PARAMS which is included in appendix I. A brief 
summary of these graph characteristics is presented below, while a more 
thorough discussion is reserved for article 3.2. 

PARAM-1 is a graph of alpha versus beta, with magnitude held at a 
fixed value while drawing constant omega curves. Interchanging alpha and 
beta produces absolutely no alteration of curves, other than a ninety de- 
gree re-orientation. PARAM-2 is a graph of alpha versus beta again, with 
omega held at a fixed value this time, while constant magnitude curves are 
drawn. As in PARAM-1, no difference other than a ninety degree re-orienta- 
tion is obtained when the axis variables are interchanged. PARAM-3 places 
alpha on the ordinate, and magnitude on the abcissa. Omega is held at a 
fixed value throughout the computational process, while constant beta 
curves are generated. Now interchanging alpha and beta in PARAM-3 will 
result in an entirely different set of curves, since the sensitivity of 
these two variables is normally different. PARAM-4 is a graph of alpha 
on the ordinate versus omega on the abcissa. Beta is held constant while 
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magnitude curves are generated. Again, different curves will be obtained 
when alpha and beta are interchanged, since holding one characteristic 
equation parameter constant will produce an entirely different effect if 
another parameter is held at a fixed value instead. PARAM-5 generates 
constant omega curves with beta held at a fixed value, and graphs these 
curves against magnitude on the abcissa, and alpha on the ordinate. As 
with PARAM-4, the PARAM-5 curves are quite different when the roles of 
alpha and beta are interchanged. PARAM-6 is a graph of alpha on the ord- 
inate and omega on the abcissa. Magnitude is held at a fixed value through- 
out the calculations while constant beta curves are generated on the alpha- 
omega plane. This is the same arrangement as PARAM-3 with the roles of 
magnitude and omega being reversed. Interchanging the alpha and the beta 
in PARAM-6 will also cause a sharp alteration in the curve configurations. 
The familiar Bode plot is herein reproduced as PARAM-7. Needing no re- 
view, it is simply a logarithmic plot of magnitude on the ordinate, and 
omega along the abcissa. With beta held fixed at some constant value, 
constant alpha curves are generated. Conversely, when alpha is held 

fixed at some constant value, the beta curves are generated and plotted 

on the magnitude - omega plane. One problem exists with the logarithmic 
plot, however. In order to draw a useful logarithmic grid, it was neces- 
sary to shift the origin to the lower right hand corner of the graph paper. 
In this manner, the decade span could cover a total of 15 inches rather 
than a maximum of nine. The net effect has been to require the abcissa 

to be plotted along the right hand edge of the graph while the ordinate 

is plotted along the bottom where the abcissa otherwise would be, thereby 
coinciding with the X-axis. Thus the ordinate originates in the bottom 
right hand corner and increases toward the lower left hand corner, or 


original origin. The same process is required for PARAM-4 and PARAM: 6 
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whenever specifying the logarithmic plots. Detailed coverage of the 
available options will be found in article 3.2 which is devoted to 
program PARAMS. 

The programs PARAM-D and PARAMS utilize only four variables, al- 
though their adaptation to others has been made. The method of solution, 
however, has been with two variables along the axes, while the remaining 
two have been employed as constant curves, and in the case of program 
PARAMS one variable has been ehld fixed at a constant value. In program 
PARAM-D, the axis variables have been selected as alpha and beta only; 
no attempt has been made to develop any other projection whereby other 
variables might replace alpha and/or beta as axis coordinates. Program 
PARAMS has attempted to focus attention on several two dimensional pro- 
jections in the hopes that peculiarities not evident in one planar 
presentation might become conspicuous in another mode. Due to the ex- 
tensive additional projections with the incorporation of one more vari- 
able, it was decided to refrain from including phase ( © ) as a para- 
meter in all but the Bode diagram, PARAM-7, where its utility and value 
have already been demonstrated. 

The following article will develop the basic capabilities of program 


PARAMS with regards to its application in frequency response studies. 


3.2 Digital Computer Program PARAMS. 

The preceding article and earlier example in section 2 presented in 
some detail the capabilities of program PARAM-D for stability studies. 
Program PARAMS has been developed to study some pertinent factors in fre- 
quency response which might prove useful in systems and controls design. 

For each of the graphical representations entitled PARAM-1 through 


PARAM-7, a separate computer program was developed in FORTRAN II language 
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for use on the Computer Data Corporation (C.D.C.) Model 1604 digital 
computer [11]. Each of these programs was compiled in a similar fashion, 
and employed the algebraic methods derived by Doctors G. J. Thaler and 
A. G. Thompson [12]. As previously discussed in this section, each graphi- 
_cal configuration has been developed in order that the controls design 
engineer might gain additional tools to facilitate studies in analysis 
and synthesis. For the study of one particular problem, it is often 
necessary to analyze the systems behavior from wany different angles, and 
then select a smaller number of characteristics for further intensive 
study. Unfortunately, every real life problem involves a different set 
of constraints, and for that reason there can be no set pattern for suc- 
cessful design application in all cases. Many situations reduce to trial 
and error techniques on preliminary studies in order to obtain the best 
compromise when a number of possible solutions exist. The ability to ac- 
complish this objective using many separate computer programs is indeed 
possible but cumbersome. Instead of this method, one consolidated com- 
puter program titled PARAMS has been written in FORTRAN IV language for 
use on the International Business Machine (IBM) model 360 computer. The 
entire computer program, together with explanatory flowgraphs and illustra- 
tive printouts, has been included as Appendix I. Details of its use, to- 
gether with the advantages and disadvantages, will now be presented. 
Versatility is indeed a virtue, provided it does not gain domination 
over simplicity. While considerable apprehension might be the first re- 
action upon reviewing Appendix I, a systemized pattern will become evident 
upon further inspection. With this program, complexity does not enter into 
the computational procedures, but rather a considerable amount of detail 
has been included in the printout format for the sake of clarity when re- 


viewing the material contained within the graphs. At most, the computational 
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process is intricate in the initializing stage wherein the various axis 
maximum and minimum limits, together with the curve and fixed constant 
values are redesignated for the ensuing common computational and graphi- 
cal construction sections. This can be verified by noting that the numera- 
tor and denominator coefficient calculation sections are common to PARAM-1 
through PARAM-6. It will further be seen that the gain computation for 
PARAM-7 employs different formulae but precisely the same logic, and for 
that reason is incorporated in the same numerator and denominator coef- 
ficient calculation sections as used for PARAM-1 through PARAM-6. Only the 
phase calculations for PARAM-7 require a different process for generation 
of the coefficients. For this reason, this computational section is by- 
passed for all other processes. 

The logical commencement point for design studies is to define a 
range of parameter values which guarantees all roots to lie in the left 
half plane, for clearly the question of stability must remain foremost in 
the designer's mind. Program PARAM-D has been devised to meet this demand, 
and its capabilities to provide a stability perspective have already been 
delineated. In determining the various possible frequency response 
characteristics for a given system, several graphical representations 
will be invariably studied. Most likely, however, program PARAM-7 will 
be employed in the final analysis as an overall check on system specifica- 
tion criteria. The remaining programs, PARAM-1 through PARAM-6 are most 
informative in the initial and intermediate stages of design study. It 
would therefore seem that the inclusion of all frequency response para- 
meter plane programs into one package would be most useful. Different 
design studies will normally emphasize certain performance criteria more 
than others. Thus seldom would the same sequence of graphical presenta- 


tions be called for in two consecutive studies. In fact, PARAM-7 might be 
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omitted for the initial coarse adjustments, but yet included in the final 
stages of refined improvements. In this manner each system study can be 
tailored to meet certain specifications and suit specific needs. Also, 
certain systems will be more complex than others. Therefore, more than 
one mathematical model might be considered for simultaneous inspection. 
The program is constructed to handle any number of such studies at one 
time. Each of these normally requires a complete set of data cards, and 

a brief explanation of their usage will follow. Additionally, an abbrevi- 
ated description section appears in the program preface, found in Appendix 
I. The total number of these data card sets must be indicated on the 
first card following the main program. This desired number of runs is 
designated by the variable NRUN. The program simply processes the re- 
maining data cards by commencing a "DO LOOP" with the variable IRUN going 
from 1 to NRUN. The first card of each data set indicates which graphi- 
cal presentations (PARAM-1 through PARAM-7) are desired. Again, there is 
no limit to the number of these. The program simply commences another 

"DO LOOP" at statement number 26, with IPARAM going from 1 to KPARAM. 

This cycles each of the graphical constructions in the order called for. 
In this manner, a great multiplicity of variations exist for the simultan- 
eous study of different preoblems. At this point it is pertinent to in- 
dicate a technique by which additional finesse may be utilized. Although 
most studies involve extensive rearrangement of a data set from one system 
to another, some require only a slight modification. In these latter cases, 
it is not necessary to submit an entire additional set of cards for each 
study. By first specifying the desired graphical plots of one system, it 
is possible to consecutively list the desired sequence for another system. 
Then a simple insertion of the desired data modification immediately follow- 


ing statement number 26 will alter the original data in the specified manner. 
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As an example, suppose one desired to first construct PARAM-5 with alpha 
on the ordinate, and then with beta on the ordinate. Simply list KPARAM 
as 2, and then indicate 5 twice for two consecutive plots of PARAM-5. 
Insert the following statement on an additional card behind statement 26: 
"IF (IPARAM. GE. 2) MODORD = 2''. Thus the required change in the input data 
is effected at the desired time by one card rather than a complete set of 
data cards. 

Although no standard names have been assigned to these different pro- 
grams the title of each graphical representation (PARAM-1 through PARAM-7) 
has been retained to eliminate confusion in terminology for those persons 
previously exposed to the partial programs which are already in existence. 
The result has been a rather complex programming sequence for printout of 
graphical information and preliminary initializing computations, but as 
has already been stated, the basic manipulations are common to all pro- 
grams, and are centrally located in the program from statement number 350 
through statement 426. A cursory inspection of these 76 statement numbers 
as well as the remainder of the program will reveal that its construction 
is modular in form, and each section is a separate entity unto itself with 
adequate illustrative sub-titles. The overall program has been so ar- 
ranged that many sections can be removed if it is desired to employ their 
use elsewhere. Conversely, additional extensions are easily facilitated 
by this manner of program fabrication. 

Careful consideration has been used in this program design to retain 
flexibility for use on multiple problem types. Although the program is 
slightly larger in size than any of its component programs which were 
previously mentioned, it is smaller than the combination of any two of 
them. While intensive effort to streamline programs is not of particular 


significance, the repetition of computer processing is needless and 
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inefficient. Also, the preparation of data cards for a large number of 
programs serves only to increase the probability of error with consequent- 
ial waste in time, expense and labor. More important is the fact that 
computer programming is still a dynamic field which often necessitates 
program modifications to meet the requirements of more sophisticated equip- 
ment. Thus one might be motivated to use a single computer program in 
lieu of eight separate ones, and to this extent consolidation is warranted. 
To ensure adequate clarity of results, an elaborate printout of in- 
formation has been established. This information can be basically cata- 
gorized into three groups, each of which is controlled by a single input 
data variable IWRITE. The first reproduces all input data for verifica- 
tion, and is effected whenever the value of IWRITE is any non-negative 
number. The second category gives detailed information about how each 
graph will be constructed. Included in this section are the graph titles, 
the range of values for the ordinate and abcissa variables, the curves to 
be generated, and the parameter that is to be held constant throughout the 
computational process of that particular program. This explanatory sec- 
tion will be printed for all positive values of IWRITE. Finally, the com- 
puted results for each curve are printed out at pre-selected intervals. 
The value assigned to IWRITE determines at what intervals the computation- 
al results will be printed, and the range of permissible values is from 
one to 900. Thus an assigned value of one would result in the reproduction 
of every calculated curve value, whereas a value of ten would give a print- 
ed value for every tenth calculated point. Approximately one third of the 
program is devoted to printing out the above information for feedback pur- 
poses. If desired, these sections can be partially or entirely bypassed 
by the appropriate selection of the input data variable IWRITE. Any nega- 


tive value assigned to this variable will provide complete calculation, 
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but without any printout whatever. 

It has already been pointed out that if the roles of alpha and beta 
are allowed to be interchanged, a much smaller number of graphs are re- 
quired. In order to accomplish this objective, some criterion must exist 
for the choice of which variable (alpha or beta) will be plotted on the 
ordinate in the case of PARAM-1 through PARAM-6, and which variable will 
be held at a fixed constant value while the other generates the desired 
curves in the case of PARAM-7. The input data variable titled MODORD has 
been designated for this purpose, and so by this selection the type of 
curves is controlled. With a value of one for MODORD, alpha curves will 
be generated in the case of PARAM-7, and alpha will be plotted on the 
ordinate for the rest of the programs. Conversely beta will be used in 
lieu of alpha for those functions when the value of MODORD is set equal 
to two. This flexibility for PARAM-1 and PARAM-2 might be unnecessary, 
because as indicated earlier, a ninety degree reorientation of the graph 
is the only effect of interchanging alpha and beta when these variables 
are on the axes. However, the present day graphical plotting facilities 
are generally constructed in such a manner that the abcissa, or X-axis, 
is more severely restricted in length than the ordinate, or Y-axis. The 
facilities available at the inception of program PARAMS were in fact such 
that the ordinate length could be approximately twice as great as the 
length of the abcissa. It was envisioned that in some cases an expansion 
of one variable might be more desirable than another. Extending the ap- 
plicability of MODORD to PARAM-1 and PARAM-2 makes this provision possible, 
which otherwise would require modification of the main program to achieve 
such results. Furthermore, it is advantageous to retain a common orienta- 
tion of graphical construction when viewing several projections simultan- 


eously. 
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The value of magnitude often occurs in several different forms, and 
thus it provides a source of confusion if the terminology is not clearly 
specified. The input data variable IMAG serves to eliminate this poten- 
tial error by distinctly categorizing magnitude into one of three forms. 

The program is constructed in such a manner that the form of magnitude as 
specified by IMAG is maintained consistently for all graph plots. Al- 
though this is strictly a matter of convenience, it was found that various 
approaches to design require magnitude to be handled in different forms; 
hence this inclusion. Thus for magnitude to be interpreted as decibels, 

a value of zero would be assigned to IMAG, and the curve values of PARAM-2 
and PARAM-4 would be listed in decibels, while the abcissa values of PARAM- 
3 and PARAM-5 would likewise be interpreted accordingly. If actual magni- 
tude values were to be used instead, IMAG would necessarily have to be 
changed to one. Assigning a value of two to IMAG would dictate all values 
to be handled as magnitude squared. In all facets of life it seems that 
most conventions are fraught with exceptions, and this program has not 
escaped that dilemma. Since all Bode plots associate magnitude with deci- 
bels, it was decided to maintain this practice; consequently the selection 
of IMAG has no effect on the form of PARAM-7 graphs. They automatically 
plot magnitude along the ordinate in the form of decibels. Conventional 
practice has also dictated that omega be logarithmically plotted along the 
abcissa for all Bode diagrams. Other frequency response graphs which also 
plot omega on one axis have not acquired any standardized form, and for that 
reason provisions have been made to select a square linear grid or a logari- 
thmic grid as deemed appropriate. The input data variable IGRAPH provides 
this option for the various possible combinations of PARAM-4, PARAM-6 

and PARAM-7. Additionally, this variable is used to prevent any graphs from 


being generated by assigning to it a value of minus two. This provision is 
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useful only for efficient use of available hardware when desired graph 
scales are initially unknown. 

In order that maximum versatility might be gained, the memory stor- 
age requirements of this program have been expanded to present day capa- 
bilities. A review of the program DIMENSION statement just preceding the 
input data read in section will show that the bulk of the storage capa- 
city is consumed by the two dimensional variables 5 through W. 

These are employed for the various optional arrangements of PARAM-7 and 
may be combined into one common storage array for all other calculations 
when employing this program in smaller computer facilities. In this 
manner, the memory core may be reduced from 500K bytes to less than 200K 
bytes. This is accomplished by inserting the following statement after 
the dimension statement: EQUIVALANCE (S,T,U,V,W). 

The additional variable, phase, has been incorporated in PARAM-7 
with several variations possible. The input data variable IPHASE is 
used to govern the relation between the magnitude and phase plots. 

No phase curves will be drawn when IPHASE is set equal to a negative num- 
ber. Setting IPHASE equal to zero will enable phase curves to be super- 
imposed upon the magnitude curves for ease in investigating stability 

by means of phase and gain margins. Assigning any positive number to 
IPHASE will provide separate phase curves. To allow for various super- 
positions, and various separate phase curve configurations as well, it 
was decided to employ an additional variable IQUAD which allowed the phase 
graphs to be plotted in the following two manners: When IQUAD has an as- 
signed value of zero, phase curves will be automatically scaled at 45 
degrees per inch, with 180 degrees on the magnitude axis. When the vari- 
able IQUAD is positive, however, the phase curves will be plotted with a 
scale equal to the value of IQUAD, starting with 0 degrees at the top of 
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the graph. In this manner, it is possible to cover any desired range of 
phase angles. 

Basic conformal mapping theory indicates that a polar plot can be 
generated from the information contained in the magnitude and phase plots. 
Particular importance might be attached to the polar plot representation 
of a system when the Nyquist criterion is preferred over the Bode stabil- 
ity analysis method. Although it is a duplicative process and therefore 
not strictly required for information content, its inclusion has been 
assumed appropriate for the sake of completeness. The basic construc- 
tion of the polar plane and polar curves are contained within the program 
statement numbers 537 and 454 under the heading "PARAM-7 Polar Plot Curves". 
Its usage is governed by the rather appropriate input data variable IPOLAR, 
and a polar plot will be drawn when its value is one. Conversely, no 
polar plot will be drawn for any other value of IPOLAR. 

An additional provision has been incorporated in program PARAMS 
for extended usage of the Bode plot. Certain advanced feedback design 
techniques call for complex systems to be studied from a cumulative stand- 
point. Such design measures require vector addition of component Bode 
diagrams and as such, the linear addition of each magnitude plot and each 
phase plot will yield erroneous results. Instead, the vector quantities 
of each design segment are summed in complex notation, and the final re- 
sult is then plotted as a single magnitude curve, and also a single phase 
curve when the latter is requested. The program subsection designated to 
accomplish this process is titled ''PARAM - 7 Superposition Section", and 
is contained by the program statement numbers 520 and 529. It is con- 
trolled by the input data variable NPAR7. where the magnitude of this 
variable dictates how many component terms are desired to be added toget- 


her. It is necessary to emphasize at this point that this section performs 
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a vector superposition of feedback terms, and not merely a scalar super- 
position. Thus whenever NPAR7 is greater than one, this section will 

add NPAR/ minus one terms vectorially to the first term before graphing 
the results. In order that these components may be properly calculated, 
provisions have been made to insert the component equation terms in the 
same coefficient format as normally submitted for the conventional cal- 
culations. These additional terms, however, are appended to the initial 
group of input data cards. Therefore NPAR7 minus one groups of coeffici- 
ents must be submitted directly behind the title cards of the rest of the 
program. 

A brief summary of the basic operational procedures which are common 
to each of the programs will now be undertaken. The schematic corollary 
in appendix I will parallel the following outline and should be consulted 
together with the main program for more detail. 

The first six statements in the "Input Data Read in Section" pro- 
vide jin proper FORTRAN language, the necessary format for most of the 
read and write statements of the program. Their use is specified in the 
"Data Card Information Section" which appears at the beginning of the 
program. The remainder of the input data read in section is devoted to 
reading all the 19, 21, or 23 data cards as outlined in the preceding in- 
formation section. The next step is to print out the input data when so 
requested for verification. This procedure is accomplished by the section 
bounded by the statement numbers seven through 23, inclusive. Following 
the input data verification section, the program initializes the various 
parameters and computes the necessary controlling variables from such in- 
put data as axes maximum and minimum values. Within this section, con- 


tained by statements 24 to 54, the desired value of NPARAM is selected and 
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this value guides the rest of the program sequence by assigning the proper 


axis variables and their limits, the correct number and type of curves to 


be generated etc. 


When the logarithmic plot is requested for PARAM - 4 


or PARAM - 6, and whenever the Bode plot is to be generated, the section 


titled "Log Grid Graphical Construction" is executed. 


Extending from 


Statement number 55 to 58, this sub-section utilizes the information 


generated in the initializing section 
and origin shift if applicable. Then 
of decades by rounding up the maximum 
of ten, and rounding down the minimum 
ten. Utilizing these factors and the 


constructed in such a manner that the 


to determine the ordinate scale 

it calculates the required number 
limit of omega to the nearest power 
omega value to its nearest power of 
axes lengths, a logarithmic grid is 


omega values are properly labelled 


along the right hand edge of the graph paper, while the ordinate values 
are indicated across the bottom in such a manner that their magnitudes are 
correct but the signs reversed. Thus the ordinate scale must be correct- 
ed to read increasing from right to left. Although this arrangement is 
less than optimum, it has proved to be much better than some alternatives 
which were previously attempted. 

The most prolific section of the program ranges from statement number 
59 to 200 under the sub heading "Graphical Information Printout Section". 
Its utility lies in the desire to establish positive correlation between 
printed data and graphical presentations when several similar graphical 
analyses are computed simultaneously. Following the extensive information 
printout, the main computational process commences by initiating calcula- 
tions for the first curve with statement number 202. Next, the coefficient 
preceding the quadratic formula square root is assigned by the variable 
ZSIGN. 


Following this designation, the abcissa is incremented 900 times 


from the minimum specified value to the maximum value. Additionally the 
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curve value is selected, based upon the value of the variable in statement 
number 202. Then the appropriate parameter which is to be held at a 

fixed constant value is assigned its pre-selected value. The ensuing 
calculations commencing with statement number 350 are in conformance with 
the theory presented in section 2. The magnitude protion of the Bode 
program requires the coefficients designated SUMN for the numerator and 
SUMD for the denominator, while the remaining programs utilize the 
numerator variables SNB2, SNB1, and SNBO in addition to the denominator 
variables SDB2, SDB1, and SDBO. Such requirements are generated under 

the sub headings "Numerator Coefficients Calculations" and 'Denominator 
Coefficient Calculations". Following these sections, the necessary phase 
coefficients are computed from statement number 386 through 400, and the 
additional Bode curve ordinate value calculations are continued through 
statement number 410. The section titled "Quadratic Formula Solution of 
Dependent Variable, PARAM-1 through PARAM-6"" is self explanatory ana 
simply computes the ordinate value commensurate with the assigned fixed 
value of the parameter which is kept constant, together with the specified 
value of the curve which is being generated, and the corresponding point 
on the abcissa. The computed results are then printed in accordance with 
the provisions contained between statement number 426 and 448. Having 
computed this ordinate value, the abcissa which is the independent vari- 
able is again incremented and a corresponding ordinate, or dependent vari- 
able point is calculated. After reaching the maximum permissible abcissa 
value, the process is again repeated for the alternate solution of the 
quadratic formula. When both possible ranges of values have been scanned, 
the next curve is specified, and the entire procedure repeated. Upon ter- 
minating the curve requirements, the program then reverts to statement 


number 26 for the next cycle in the run. When all requested graphical 
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construction cycles have been accomplished, the next run is processed 
in like manner. The overall program is terminated when all desired runs 
have been computed and graphed. 

The graphical plotting procedures will be discussed separately be- 
cause their intricate operations tend to destroy the continuity of the 
preceding discussion on the overall program operation. Basically, the 
graphical plotting subroutine is executed whenever the following condi- 
tions are engaged: 1) upon initially setting up the logarithmic plot 
when specified in accordance with the above dictates; 2) when two or more 
computed points fall within the prescribed graph range and the next point 
lies outside the ordinate graph limits; 3) when the end of each scan of 
the abcissa has been reached - once for each value of the quadratic equa- 
tion square root coefficient; 4) at the end of each set of curves for one 
graph, as accomplished by the graphical termination section found at the 
end of the program beginning with statement number 999. 

The versatile arrangement necessitates several provisions for the 
many possible combinations of the Bode plot. Since the quadratic formula 
is not employed in these calculations, only one incrementing scan of the 
abcissa is performed. Additionally, the computed magnitude values have 
been programmed in such a manner that when the ordinate graph limits are 
exceeded, the external points are fictitiously redefined to lie on the 
appropriate ordinate graph limits. This procedure was necessitated by the 
additional phase calculations which remained within the graph limits 
throughout the entire range of the prescribed omega values. Finally, a 
separate provision had to be made for the construction of the phase curves, 
since an additional logarithmic graph was required when phase curves were 
plotted separately. This procedure has been accordingly programmed in the 
section titled "PARAM - 7 Phase Curves" found between statement numbers 529 
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and 537. The section titled ''PARAM - 7 Polar Plot Curves” has already 
been discussed and is a separate entity unto itself, utilizing the pre- 
computed values of magnitude and phase, and plotting same in polar coordi- 
nates. This process is accomplished between statement numbers 537 and 999 
in the program. 

The preceding development has summarily shown that this program con- 
tains broad applications, and is well suited for additional modifica- 
tions to meet particular needs. Blanket provisions for another graphical 
configuration has been purposely included in the program sequence labelled 
NPARAM- 8, and should future expansion be deemed appropriate, the required 
modification will merely entail filling the appropriate spaces with the 
desired initializing variables and applicable computational equations. 

To facilitate such extensions and also to enable maximum comprehension 
of the program operations, two logic flow diagrams have been included in 
Appendix I. The first is a generalized sequentail listing of logical 
operations which are applicable tc all graphical configurations. The 
second is a detailed flow chart for PARAM - 6, so chosen because it in- 
cludes most of the alternatives presented in this discussion. 

A final comment is appropriate concerning the application of computer 
program PARAMS. As indicated earlier in this article, considerable atten- 
tion has been given to versatility. The associated price of computational 
time and memory storage has consequently been paid. The procedures to be 
used in reducing memory storage for PARAM-1 through PARAM-6 has already 
been indicated. Compilation time is roughly 1.5 minutes. An additional 
30 seconds is consumed for each printout of preliminary graphical informa- 
tion. The actual calculation of 16 curves for each graph requires approxi- 
mately 3 minutes. Thus a considerable amount of computer time is consumed 


when a complete set of curves is requested for all possible graphical 
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representations. Various methods exist to reduce this time requirement. 
By submitted an object deck, compilation time is deleted. This imposes 
the restriction of submitting a complete data card set for each run, 
however. Eliminating the input data printout, and the preliminary graphi- 
cal information printouts will also reduce the time requirements some- 
what. This, however, is not recommended because the elaborate informa- 
tion presentation provides an exceedingly valuable check on the validity 
of the graphical results. Finally, the most time is consumed by a highly 
sensitive incrementation of the abcissa variable. It will be recalled 
that the range of the abcissa variable is divided into 900 evenly spaced 
steps. For each of these abcissa values, a corresponding ordinate value 
is computed. Although many situations lend themselves to a larger in- 
crementation, it was decided to maintain the fine step sizes for highly 
sensitive system curves. The maximum sensitivity of 900 points per graph 
line is the restriction placed upon the program by the CALCOMP plotter in 
use at the Naval Postgraduate School computer facility. Reduction of this 
incrementation sensitivity is therefore another means by which computation- 
al time can be further reduced. The associated loss of information con- 
tent must be considered in such cases. As will be pointed out in the 

next section, not all PARAMS frequency response curves would be required 
or desired in every design study. The judicious choice of the most valu- 
able plots for a given set of constraints will normally eliminate approxi- 
mately half of the program possibilities. Only by experience, insight and 
intelligent design procedures can a good trade-off be achieved between the 


tools available and the data required. 
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4, Tachometer and Acceleration Feedback Compensation Studies. 

The mathematical developments of section 2 have provided the basis 
for constructing the computer programs which were discussed in Section 3. 
This section will apply these methods to a system study employing tacho- 
meter and acceleration feedback as compensation techniques. The frequency 
response parameter curves will be analyzed according to the respective 
utility of each as a design tool. Meticulous detail is employed in the 
initial discussion of these curves. This is to provide a clear under- 
standing of the limitations and advantages of each respective graphical 


representation. 


4.1 Third order 

Although an entire library of control engineering studies has been 
created, there has been no attempt to distinguish a set of example prob- 
lems as superlative for generalized case studies. Therefore a simple 
third order system has been considered for an initial application of the 
parameter plane frequency response curves. The illustrative example which 
was introduced in section 2 will now be explored in depth. It will be re- 
called that Figure 2-3 presented the schematic diagram of the third order 
system with compensation included. The stability curves for the system 
where mapped into the parameter plane shown in Figure 2-5. In that pre- 
sentation alpha, the abcissa variable, represented acceleration feedback 
( us ). Similarly beta, the ordinate variable, indicated tachometer 
feedback ( ee ). Reference to the characteristic equation of the sys- 
tem shown in Figure 2-3 shows that the steady state velocity error is 


aeRO, +e Kale _ 500 + 100, 000¢ 


K 
si K 100,000 
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Inspection of Ky shows that the smallest allowable value of beta 
will minimize the steady state velocity error. It is presumed that this 
is one of the predominant performance criteria. Furthermore, one can 
arbitrarily stipulate the dominant complex roots to have a damping ratio 
( 5 ) of 0.5. With these system stability requirements, one can see 
that optimum performance will occur at point A in Figure 2-5. Thus the 
optimal value of beta is 0.036, while the optimal value of alpha is 
0.00047. 

A tacit assumption up to this point has been that frequency perforn- 
ance constraints are of minor importance. Most design studies do not 
conform to this type of ommission; however stability is of interest in 
all systems. It has therefore been assumed that the two stability per- 
formance criteria above are desirable but not mandatory. This permits 
frequency response specifications to be considered with some initial 
reference point. Figure 4-1 shows an approximate correlation between 
frequency response, transient performance, and root locations in the com- 
plex frequency plane. Additionally the separate root loci with alpha and 
beta varied in Figures 4-2 and 4-3 respectively show the approximate vari- 
ation of root locations as a function of one parameter only. The corres- 
ponding effect in the frequency domain is represented by the Bode plots 
in Figures 4-4 through 4-7. The graphical relationships presented in 
Figures 4-1 through 4-7 are certainly not new to the design engineer. 
Rather, they are the fundamental tools used in most system studies today. 
As such, they are included primarily to establish a firm foundation upon 
which parameter plane frequency response curves can be built. In this 
manner, a better understanding can be gained concerning the utility of 


the latter curves which are unfamiliar to most engineers. 


94 


PARAM - 1 

Complete analysis of the information contained in PARAM-1 plots 
will now be endeavored. Figure 4-8 shows the relationships between alpha, 
beta and omega which exist for a constant magnitude of 12 decibels. It 
is clear from the legend at the bottom of the graph that alpha is the 
abcissa variable, while beta values are indicated on the ordinate. Each 
curve in the alpha - beta plane represents a fixed magnitude of 12 deci- 
bels, but a different value of omega. Thus the curve of W = 10 rad/sec 
barely appears in the lower right hand corner. As omega increases, the 
curves can be seen to shift toward higher values of beta and lower 
values of alpha. It can also be seen that the area enclosed by each 
curve decreases as omega increases to W = 70 rad/sec. The meaning of 
each of these 14 curves can be explained as follows: For the third order 
system in question, any point on a specified omega curve depicts the 
locus of alpha and beta values which will guarantee a magnitude of 12 
decibels at that specific frequency. This by itself does not portray 
much information, because it says nothing about the area enclosed by the 
curves, nor the area outside the locus of these curves. In most situations, 
it is possible to glean considerably more information out of these curves. 
Consider for a moment that value of alpha and beta at which two omega 
curves intersect. Point B of Figure 4-8 shows one intersection of the Y= 25 
rad/sec and W = 30 rad/sec curves. For this value of alpha and beta, a 
magnitude of 12 decibels is guaranteed to exist at both values of omega 
indicated. This observation is important, for it means that the magni- 
tude is either greater than 12 decibels, or less than that value for the 
intermediate range of omega. Furthermore, it is possible to assume that 


the magnitude cannot be both greater than and less than 12 decibels for 


35 


these intervening values of omega. The reason for this is simply that 

the magnitude curve is continuous, and must pass through 12 decibels in 
order to change from a high value to a lower value, or vica versa. If a 
sufficiently large number of omega curves are drawn in the region of 
interest, it can be seen that, at Point B, only two omega curves do actual- 
ly intersect. Therefore, this range of omega values represents the 12 
decibel bandwidth of a system which satisfies these values of alpha and 
beta. With this insight, it can be concluded that a larger bandwidth is 
obtained for a point slightly below point B. Thus at a value of alpha 

and beta indicated by point C, interpolation of the omega curves shows 

an approximate bandwidth from W = 24 to W = 31. Conversely, proced- 
ing in the opposite direction, the bandwidth is reduced to the point D 

at which only one frequency curve satisfies the specified value of alpha 
and beta. No other frequency curves will intersect this point, for it 
lies on the tangent locus of all frequency curves evaluated at 12 deci- 
bels magnitude. Thus this tangent locus defines all the values of alpha 
and beta which guarantee a peak resonant magnitude of 12 decibels. The 
value of the resonant peak frequency is not specified by the tangent locus. 
The reason for this is because omega is an implicit function of the tan- 
gent locus equation. Article 2.5 discussed various methods to generate 

the required non-parametric form of this M p ( w) equation, and the 
limitations encountered in each of those approaches. For those reasons, 
the tangent locus curve was not constructed. Instead, a sufficiently 

large number of omega curves were constructed to show the approximate loca- 
tion of the tangent locus. Although the resonant peak frequency is not 
directly obtainable from the tangent locus, it is clearly seen that such 


information is readily found by the intersection of a particular omega 
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curve with this tangent line. 

The discussion thus far has presumed that all values inside the 
tangent locus represent peak resonant points greater than 12 decibels. 
This is intuitively obvious, but no proof has been made of this fact. 
Figures 4-9 through 4-12 have been included to verify this presumption, 
and also to illustrate additional phenomena. Figure 4-9 clearly shows 
that when the magnitude is reduced to 6 decibels, a greater area is en- 
closed by each of the omega curves. This results in a tangent line rep- 
resenting the locus of peak resonant magnitude equal to 6 decibels. 
Comparison of this tangent locus with the one in Figure 4-8 will reveal 
that the 12 decibel tangent locus falls within the area bounded by the 
omega curves of Figure 4-9. The points B and D of Figure 4-8 have been 
reproduced on the constant 6 decibel graph for comparison. It can be 
noted that point B represents a 12 decibel bandwidth of 5 rad/sec and a 
6 decibel bandwidth of approximately 13 rad/sec. Additionally, point D 
depicts a resonant peak magnitude of 12 decibels at 28 rad/sec with a 6 
decibel bandwidth of 13 rad/sec. Similar information can be obtained 
for other values of alpha and beta throughout the parameter plane. By 
similar analogy, it is possible to determine what range of alpha - beta 
values will guarantee a peak resonance not to exceed a pre-selected value. 
Thus all values between the M Pr (w) = 12 and the M p (w) = 6 curves 
will guarantee all magnitudes to be less than 12 decibels, with a 6 
decibel bandwidth of variable size depending on the choice of alpha and 
beta. A closer inspection of Figure 4-9 will reveal that the maximum 6 
decibel bandwidth attainable is approximately 17 rad/sec for the situation 
just described. This maximum bandwidth occurs approximately alpha equal 
to 0.0001 and beta equal to 0.018. At this point the resonant peak magni- 


tude is slightly less than 12 decibels, with a resonant angular frequency 
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of approximately 40 rad/sec. The cut-off frequencies occur at W = 32 


rad/sec and W = 49 rad/sec. All of the above information has been 
derived from the evaluation of two PARAM-1 graphical plots. It will be 
recalled that the PARAM-7 curves of Figures 4-4 through 4-7 merely por- 
trayed Bode curves as a function of one parameter only. These PARAM-1 
curves denote a portion of many Bode curves as a function of two system 
parameters. The value of these two is restricted to only two magnitudes, 
but certainly the analysis of any two such graphs provides a certain in- 
sight as to the system behavior for other magnitudes. Interpolation can 
be used for the intermediate magnitude range. Such a procedure is not 
recommended without some caution for extensions of magnitude studies out- 
side the graphed limits. The consideration of Figures 4-10 through 4-12 
will illustrate this point more clearly. For a constant magnitude of 3 
decibels, Figure 4-10 shows a dramatic increase in the applicable range 

of alpha - beta values. The graph scales of all PARAM-1 plots were kept 
the same simply to illustrate the large increase in area encompassed by 
each constant frequency curve as the value of magnitude is decreased. 

As might be anticipated, a much larger 3 decibel bandwidth exists for 

the points, B, C and D shown on Figure 4-8. The bandwidth of these three 
points do not vary greatly, however. Each is approximately 20 rad/sec. 
This illustrates the sensitivity of beta, the tachometer feedback in the 
area under discussion. For a given value of alpha equal to 0.008, a 

slight change in beta greatly affects the size of the peak reasonant 
magnitude. Conversely, an equal change in alpha reflects a vast difference 
in omega, the location of this peak resonant magnitude, plus a slight 
change in its magnitude value. Although the sensitivity of the alpha vari- 
able might not seem apparent by looking at the curves, it must be remembered 


that there is a different graph scale for each of the axis variables. 
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Figures 4-11 and 4-12 show most clearly why the analysis of merely 
two PARAM-1 frequency response representations might be misleading. 
Consider point B of Figure 4-8 which is reproduced on the zero decibel 
magnitude graph of Figure 4-11. The constant omega curves of this plot 
show that no two curves intersect at point B. Clearly, it is seen that 
all values of omega less than 40 rad/sec encompass this choice of alpha 
and beta. Therefore a cut-off frequency of approximately 39 rad/sec is 
obtained for the zero decibel bandwidth of that system. Similarly, 
from Figure 4-12, the -3 decibel bandwidth is about 43 rad/sec. Using 
only these latter two PARAM-1 plots, one cannot discern any resonant peak 
for the system frequency response associated with these values of alpha 
and beta. It is only through the apriori information portrayed in any 
one of the Figures 4-8 through 4-10 by which we are able to make such a 
statement that a resonant magnitude actually exists. Only by the combined 
interpolation of all the Figures 4-8 through 4-10, can we state with a 
fair degree of confidence that the value of this resonant magnitude is 
approximately 15 decibels. This of course could be checked by construct- 
ing another PARAM-1 curve, or perhaps a Bode plot. Neither of these are 
necessary for initial coarse approximations. It will soon be shown that 
the combined usage of PARAM-1 and PARAM-2 provides a much better method. 
The reason for the inclusion of so many PARAM-1 curves was simply to il- 
lustrate how a valuable insight can be obtained for system frequency re- 
sponse studies when considering two parameter variations simultaneously. 
Within the limitations just indicated, it is possible to determine how 
two parameter variations affect bandwidth and peak resonance from a care- 
ful inspection of several PARAM-1 graphs. 

Before proceding to PARAM-2 graphical configurations, it will be 


advantageous to review Figures 4-8 through 4-12 in a three dimensional 
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concept. Consider for a moment a curved conical surface similar in 

shape to a bull's horn. Also, let this curved cone be elliptical in 
shape, rather than round. Next assume that the base of this surface is 
in the lower right hand corner, nearest to the observer. The smaller end 
of this surface would therefore be in the upper left hand corner of the 
observation sector, at a certain distance away from the observer. Now 
superimpose the omega axis on the centerline of this conical surface in 
such a manner that the low values of omega are at the base of this sur- 
face, or nearest to the observation point. The higher values of omega 
would consequently lie further away from this observation point, in- 
creasing in value toward the curved conical tip. With this three dimen- 
sional surface concept in mind, look at Figure 4-8. The curved omega 
axis originates in the lower right hand corner and extends through the 
alpha - beta plane, curving upward and slightly to the right. Figure 4-9 
depicts this curved conical surface as larger in size, (i.e., closer to 
the observer), with the base of this surface shifted slightly to the left, 
thereby bringing it more into the line of sight. Bringing this curved 
conical surface still closer to the observation saint makes it appear 
larger, as seen in Figure 4-10. The zero decibel PARAM-1 graph, Figure 
4-11, portrays an even larger shape with the base shifted even more to 
the left. Finally Figure 4-12 can be regarded observing this curved sur- 
face with its base very close, similar in appearance to passing through 

a curved tunnel under a river. 

The preceding illustrative example can be considered as elementary, 
but is fundamental to the understanding of the various graphical rep- 
resentations to follow. No attempt will be made to correlate the il- 
lustration to the PARAM-7 curves of Figures 4-4 through 4-7, since those 


are essentially single parameter variations of the series of two parameter 
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surfaces just described. Better use of the illustrative example is made 


by correlating PARAM-1 with PARAM-2. 


PARAM-2. 

Consider now a cross sectional cut of the curved conical surface 
described above. PARAM-2 effectively maps constant magnitude curves onto 
the alpha - beta parameter plane for a fixed value of angular frequency. 
In this manner, it is possible to show how a particular omega value, 
such as W = 10 rad/sec will vary as a function of magnitude in the 
alpha - beta plane. Consider again the W = 10 rad/sec curve in Figures 
4-8 through 4-11. The associated PARAM-2 plot is shown in Figure 4-13. 
Since the W = 10 rad/sec curve represented a cross section of the curved 
conical surface near its base, it is possible to observe how this base 
is rotated toward the left as magnitude is decreased. Further inspec- 
tion of Figure 4-13 shows how the curve increases in size with a reduc- 
tion of magnitude. The -3 decibel, W = 10 rad/sec curve is in fact 
shifted completely off the graph of Figure 4-13. This effectively shows 
that no alpha - beta value within the prescribed range will have a magni- 
tude of -3 decibels at an angular frequency of 10 rad/sec. No information 
is conveyed about any other angular frequency in Figure 4-13. To obtain 
such information would necessitate a different PARAM-2 plot, or a compari- 
son with one or more PARAM-1 curves. This does not preclude the correla- 
tion of other frequency response curves. They can be consulted for pert- 
inent relationships, and will be discussed later. Figure 4-14 is a graphi- 
cal representation of constant magnitude curves for a fixed angular fre- 
quency of 20 rad/sec. Inspection of these curves shows that they lie very 
closely to the centroid of the curved conical surface discussed in the il- 


lustrative example of this article. This can be verified by their 
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concentric shape. Due to the graphical construction technique, these 
curves do not extend to the bottom of the graph. This can be explained 
by the fact that these curves are almost vertical at the lower graph 

edge. Even though the abcissa variable is incremented 900 times during 
each curve computation, this is not sensitive enough to allow the detec- 
tion of such a steep slope. If an infinitessimally small step size were 
possible, these curve values would be detected. This limitation of the 
graphical plotting procedure must be kept in mind. The mathematical 
theory of section 2 indicated that two quadratic formula solutions were 
necessary for the complete solution of any curve generated by PARAM-1 
through PARAM-6. The computer technique discussion of section 3 dis- 
cussed how each graph curve was constructed first for one sign of the 
square root coefficient in the quadratic formula, and then for the other. 
Subsequent graphs will show how these two curve solutions do not always 
connect. The reason for this again lies in the finite incrementation of 
the abcissa variable. Returning once again to Figure 4-14, it is possible 
to note certain characteristics regarding the angular frequency of 20 rad/ 
sec as a function of the two parameters alpha and beta. The intersection 
of the peak values of each magnitude curve defines a vertical line in the 
parameter plane. Therefore it is possible to conclude that variations in 
magnitude for an angular frequency of 20 rad/sec are highly sensitive to 
beta values, and quite insensitive to changes in alpha. This says nothing 
about resonance magnitudes. Rather, each curve guarantees a particular 
magnitude value for a given angular frequency. Also, each curve defines 
an alpha - beta region within which the magnitude will always be greater 
than the boundary value for that value of omega. As a further enalysis of 


the system response for point B of Figure 4-8, it can be seen that such a 
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system will have a magnitude of 6 decibels at a frequency of 20 rad/sec. 
Figures 4-15 and 4-16 show considerably more information than the previous 
two PARAM-2 graphs, because the center of these eliptical curves is indi- 
cated. This center corresponds to the centerline of the three dimensional 
curved conical surface which has been repeatedly referenced up to this 
point. It should be observed that the center of these ellipses in Figures 
4-15 and 4-16 denote the maximum magnitudes attainable at their respec- 
tive angular frequencies. Thus by analogy, the centerline of the three 
dimensional surfaces in Figures 4-8 through 4-16 serves a dual purpose. 
Not only does it specify a particular angular frequency, but it denotes 
the locus of the highest peak resonant magnitudes for any alpha - beta 
combination of the system! This is an important conclusion, because it 
indicates whether a desired peak resonant magnitude value can ever be 
reached for the system under consideration. The comparison of PARAM-1 

and PARAM-2 graphs further shows how the conical surfaces expand about 
this centerline with decreasing magnitude. From these, it is possible 

to visualize the locus of constant magnitude conical surfaces, each 
surrounding the centerline just discussed. It is therefore seen that a 
reduction of tachometer gain coupled with an increase in acceleration 

gain will shift the resonance peak to a higher frequency. Point B of 
Figure 4-8 has also been reproduced on Figures 4-15 and 4-16. These 

show that a system with the indicated alpha - beta values will have a 
nagnitude of 12 decibels at W= 30 rad/sec, and -1 decibels at an angular 
frequency of 40 rad/sec. At this point, enough information has been ob- 
tained to obtain a fairly accurate Bode plot for the parameter values given 
at point B. Further comparison of the 15 decibel magnitude curve in 
Figures 4-14 through 4-16 show that point B does in fact have a magnitude 


equal to 15 decibels. This confirms the assumption of this value based on 
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interpolation of the PARAM-1 curves alone. Since enough information has 
been depicted in the PARAM-1 and PARAM-2 curves to portray a Bode plot 
of point B, it is obvious that similar Bode plots could be obtained for 
any other desired combination of alpha-beta points. This is exceedingly 
valuable, because it allows the design engineer to construct any number 
of coarse Bode plots by evaluating the data contained in Figures 4-8 
through 4-16. Therefore it is concluded that the combined utilization 
of PARAM-1 and PARAM-2 is a fundamental starting point for all two para- 
meter frequency response design studies. 

The preceding discussion concerning the use of PARAM-1 and PARAM-2 
was very comprehensive because of the following three key factors: First, 
a thorough understanding of new concepts can only be gained by a slow and 
detailed presentation of the uses and limitations of such innovations. 
Secondly, the combined application of PARAM-1 and PARAM-2 plots provides 
the broadest information spectrum for frequency response design studies. 
Finally, a thorough understanding of these two parameter plots is pre- 
requisite to visualizing the following frequency response projections 
which are essentially a function of one system parameter. The remainder 
of this article will cover the various roles played by PARAM-3 through 


PARAM-6. 


PARAM-3. 

Section 3 developed the major graphical presentations which were 
considered to be of most use to design engineer. In that discussion, 
it was indicated that two basically different configurations were avail- 
able in PARAM-3. Figure 3-1 showed in schematic form these differences. 
Attention will first be focused on the MODORD=1 variety as shown in 


Figure 4-17. Before discussing the information content of this graph, 
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an apparent mistake in the graphical plotting procedure is noteworthy. 
This set of curves was purposely offset to illustrate a limitation en- 
countered by careless selection of axis variable ranges. Graph scales 
are limited to one significant figure. As such the automatic scaling of 
the program PARAMS does not always produce the values which might be ex- 
pected. The requested range of alpha values was from -0.001 to 0.007. 
Since the ordinate length was specified as ten inches, the program comput- 
ed a grid scale of 0.0008 units per inch for alpha. Accordingly this 
scale was used; however, it placed the origin one inch above the minimum 
value in accordance with program plotting limitations. Thus the requested 
range of alpha values was plotted properly, even though the choice of 
ordinate values was improperly chosen. The preceding was for information 
purposes only; PARAM-3 plots will now be studied. 

Figure 4-17 portrays a set of constant beta curves on the magnitude 
- alpha plane for a given fixed angular frequency value of W = 30 rad/ 
sec. This graph can best be understood by comparing it with the PARAM-2 
plot of the same constant omega, namely Figure 4-15. Effectively, the 
PARAM-3, MODORD=1 curves take the information contained in the horizontal 
cross-sections of PARAM-2 and map this data onto the magnitude - alpha 
plane. This can be most easily visualized by taking the base of Figure 
4-15 and comparing it with the beta = 0.0 curve in Figure 4-17. With 
these two viewpoints, one readily sees another facet of PARAM-2 plots. 
Consider for a moment a topographical map of land terrain. Such maps 
represent constant altitude as a continuous curve. Each such curve is 
always a closed contour, with either peaks or valleys contained within 
the curve boundaries. PARAM-2 is such a topographical plot, with magni- 
tude equal to elevation. This concept is borne out by the magnitude- 


alpha plane of Figure 4-17. The constant beta curves of this latter 


105 


graph can be regarded as vertical slices into the topographical hill- 
side. As such, one sees which cross sections have actually been sliced 
out of the PARAM-2 magnitude contours. The complete range of PARAM-2 
alpha and beta values has been spanned by the ordinate range and constant 
beta curves in Figure 4-17. Accordingly one might conclude that a good 
indication of the system's magnitude behavior has been conveyed by the 
PARAM-3, MODORD=1 graph for W = 30 rad/sec. This is most emphatically 
not true! Analysis of Figure 4-17 indicates that the greatest magnitude 
occurs at beta = 0.0 and that the magnitude contours decrease as beta ap- 
proaches 0.05. Careful review of Figure 4-15 will disclose a serious 
fallacy in such a conclusion. Figure 4-17 shows no indication of the 
magnitude peak resonance at beta = 0.004. Without this apriori knowledge, 
a completely erroneous interpretation would have been made by analyzing 
the curves in PARAM-3, MODORD=1. This example stresses once again the 
importance of viewing system frequency response as a function of the two 
system parameters. 

Consideration will now be given to the alternate PARAM-3 curve desig- 
nated by MODORD=2. This graph presents constant alpha curves in the beta 
magnitude plane for a fixed omega. As such, these curves simply slice 
sane cross-sectional areas out of the magnitude contours of PARAM-2. 
Figure 4-18 shows such a set of curves for a constant angular frequency of 
30 rad/sec. Looking at Figure 4-15 again, all curve values lie within a 
fairly small range of alpha values. Accordingly, one would only expect 
graph curves for alpha less than 0.0021 to lie in the beta - magnitude 
plane for magnitudes greater than or equal to -3 decibels. Furthermore 


one would expect a series of magnitude peaks to occur at beta = .004 with 


the maximum amgnitude occuring on the alpha = 0.0005 curve. These factors 


0.0005 curve has maximum 


are verified in Figure 4-18 where the alpha 
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magnitude greatly in excess of 21 decibels. All other curves represent 
lower magnitude ranges as was anticipated. In this set of alpha curves, 
the maximum magnitude curve was intentionally constructed to indicate 
the large magnitude at peak resonance. This curve might have been over- 
looked, were it not for the apriori knowledge gained from Figure 4-15. 

The preceding discussion of PARAM-3 curves showed serious limita- 
tions in their usage when the PARAM-2 curve of the same constant omega is 
not utilized. When used in conjunction with PARAM-2 graphs, however, the 
PARAM-3 curves do provide some useful information. They are particularly 
useful for the erea enclosed by the largest value of magnitude in the 
PARAM-2 plots. This area tells nothing about the peak resonant magnitude, 
other than the frequency at which it occurs. By constructing one PARAM-3 
curve along the bisector of these magnitude curves, one quickly gains the 
value of peak magnitude and the magnitude sensitivity about this resonant 
point. Should the precise peak magnitude location not be so obvious on 
other PARAM-2 plots, a series of curves in the PARAM-3 graph will clearly 
indicate its location. Due to the elliptical shape of the PARAM-2 magni- 
tude curves for this system, the PARAM-3, MODORD=2 configuration would 
be better suited than the PARAM-3, MODORD=1 graph. The reason for this 
is simply that the magnitude is more sensitive to alpha than beta. There- 
fore a much smaller range of alpha values need be considered than would be 
required for beta curves. 

Up to this point, only one value of omega has been discussed. The 
PARAM-3 curves for other angular frequencies can be visualized by studying 
Figures 4-13, 4-14, and 4-16. The PARAM-2 plot of Figure 4-13 shows a 
very gentle magnitude slope along constant alpha and constant beta lines. 


Thus for constant alpha lines, one would expect the magnitude to decrease 
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slowly with an increase in beta. The magnitude would decrease more 


rapidly for larger constant alpha curve values than they would for smaller 
curve values. Conversely, for constant beta lines, it is seen that magni- 
tude will increase slowly as alpha is increased. Large constant beta 
curve values would show a slower increase in magnitude than small beta 
curve values. Figures 4-19 and 4-20 are included for verification of 
these trends. 

Inspection of Figure 4-14 shows that the PARAM-3, MODOR=1 graph will 
have a maximum magnitude at alpha = 0.0018 for all beta curves. The maxi- 
mum value of magnitude for anyvwalue of alpha will be obtained on the beta 
= 0 curve. All other beta curves will lie inside this magnitude boundary. 
The shpae of these curves will be similar to those shown on Figure 4-18. 
Since the magnitude is less sensitive at 20 rad/sec than at 30 rad/sec, the 
curves will be less pointed. The PARAM-3, MODORD-2 curves will not be of 
the shape just discussed. This is because the vertical cross-sections of 
the PARAM-2 plot show a decrease in magnitude as beta increases. All alpha 
curves will lie inside the constant alpha curve equal to 0.0018, for this 
curve represents a ridge of maximum magnitude. The constant ( = 40 Peay 
sec magnitude curves of Figure 4-16 show that the peak resonance occurs on 
the beta axis at an alpha value of approximately 0.012. The sensitivity 
of magnitude is greatest in this graph. Accordingly, the PARAM-3, MODORD=1 
curves will be much more pointed than those shown in Figure 4-17. Also 
the center of these magnitude peaks will be shifted down to the abcissa. 
The PARAM-3, MODORD=2 curves will be likewise sharper than tow shown in 
Figure 4-18, but this time the centerline of the magnitude peaks will lie 


on the beta = 0.012 line. 
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It has been shown that for this system the PARAM-3 curves are of 
limited use. Detailed correlation between PARAM-2 and PARAM-3 plots was 
made to show the inter-relationships which exist. The reason for discus- 
sing PARAM-3 graphs was to show how such information can be picked off of 
the PARAM-2 plots, thereby eliminating the requirement for another graphi- 
cal presentation. A more important reason for this detailed analysis was 
to point out the correlations applicable to this system under study. Many 
systems will undoubtedly be more complex than this. Such cases might re- 
quire the use of these curves. Therefore it is essential to obtain a firm 


understanding of these graphical advantages and limitations. 


PARAM-4. 

The graphical concepts of PARAM-4 will now be discussed for the third 
order system under consideration. Again, two basically different plots are 
obtained for this program. The PARAM-4, MODORD=1 configuration graphs 
constant magnitude curves onto the alpha - omega plane for a fixed value 
of beta. Figure 4-21 shows one such plot for a given fixed beta value of 
0.036. At first glance, these curves seem very strange. Once again, how- 
ever, recourse to the two parameter frequency response curves of PARAM-1 
will show a striking relationship. Recall for a moment the three dimen- 
sional concept that was introduced in the study of PARAM-1 graphs. Each 
three dimensional curved elliptical cone shaped surface was a function of 
alpha, beta, and omega for a given value of magnitude. As magnitude was 
decreased, the conical surface base was observed to rotate into the obser- 
vation line of sight, while the overall size expanded. The simultaneous 
rotation and expansion caused the base portion to expand more than the 
smaller end. Furthermore, recall that omega increased in value from the 


base toward the small end. The PARAM-4, MODORD=1 graph effectively slices 
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a horizontal cross section out of this PARAM-1 three dimensional sur- 

face at the prescribed value of beta. To see this more clearly, Figures 
4-8 through 4-12 will be reviewed in order. Consider first the constant 
decibel magnitude PARAM-1 graph of Figure 4-8. By cutting a perpendicu- 
lar cross section along the beta = 0.036 line, it is seen that only one 
omega curve is severed namely the W = 65 rad/sec curve. Note, however, 
that omega is a continuous elliptical surface; only a small number of 

cross sections are depicted in Figure 4-8. By interpolation, it can be 
seen that a perpendicular cross section of PARAM-1 will intersect the 

range of omega values between approximately 62 through 67. Furthermore 
note that no other omega curves have been severed in this manner. This 
intersection occurs at a range of alpha values between -0.002 and -0.004. 
Now turn to Figure 4-21 and inspect the 12 decibel magnitude curves re- 
produced there. By a similar analysis of Figures 4-9 through 4-12, it is 
possible to establish a positive correlation between the PARAM-1 and PARAM- 
4, MODORD=1 graphs. Next consider the PARAM-4, MODORD=2 plot shown in 
Figure 4-22. This configuration graphs constant magnitude curves onto 

the beta - omega plane for a fixed value of alpha. By methods analogous 

to those presented above, it is possible to slice a vertical cross section 
of the PARAM-1 graph at the prescribed value of alpha. The information 
thus obtained is seen to appear as PARAM-4, MODORD=2. The PARAM-4 curves 
are effectively topographical maps with magnitude equivalent to elevation. 
Their utility lies in the ability to correlate magnitude, frequency and one 
System parameter. They are restricted by the second system parameter which 
has been held at a constant value. As such they would only be suited for 
supplying additional frequency response information after an initial choice 
of one variable has been made. Providing a sufficient number of magnitude 


curves have been constructed on these graphs, they become in essence a 
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locus of Bode plots. When used in this manner, they become valuable 
tools. The constant alpha and beta values chosen for figures 4-21 and 
4-22 were so chosen because they coincided with the values used in 
figures 4-4 through 4-7. Thus a comparison of these PARAM-4 curves with 
the familiar Bode magnitude plots can easily be made. The PARAM-4 
curves can therefore be regarded as most useful for coarse Bode plot 
studies, since they portray magnitude-frequency information for a much 
larger number of alpha or beta values than would be feasible in a con- 
ventional Bode plot study. In this manner a design engineer would most 
likely construct encugh PARAM-1 and PARAM-2 plots to make a judicuous 
first approximation for each of the system i eaaesaae Using these, 

he could then construct both PARAM-4 plots and check his system perfor- 
mance with a fair degree of accuracy. After suitable corrections, he 
would then use PARAM-/7 for a final refined system frequency response 


performance check. 


PARAM-5 

The discussions of PARAM-3 and PARAM-4 have considered cross section- 
al cuts of the PARAM-2 and PARAM-1 planes respectively. A slightly dif- 
ferent concept is portrayed in PARAM-5 and PARAM-6. Consider first the 
PARAM-5, MODORD=1 graph shown in figure 4-23. This shows a series of 
omega curves plotted upon the alpha magnitude plane for a fixed value 
of beta = 0.025. Each constant magnitude depicts a horizontal cross 
section of the PARAM-1 graph corresponding to that constant magnitude. 
The location of the horiztonal cross-section is fixed by the constant 
beta = 0.025 for the PARAM-5 plot. This is most clearly shown by the 
following comparison. Arbitrarily choose the constant magnitude line of 


6 decibels. This automatically selectes the PARAM-1 plot shown in 


Pee 


figure 4-9. Since beta is fixed at a value of 0.025, this is the 
horizontal line of figure 4-9 to be considered. Positive correlation 
between the omega curves crossing this beta = 0.025 line in figure 4.9 
is thereby made with those omega curves crossing the magnitude line of 
6 decibels in figure 4-23. Effectively the PARAM-5, MODORD graph shows 
the series of PARAM-1 plots for any given constant beta value. Merely 
studying the different magnitude lines, one can see how a constant beta 
cross section would look on different PARAM-1 plots. The PARAM-5, 
MODORD=2 graphs show this same behavior, with the roles of alpha and 
beta interchanged. 

A comparison can also be made with a series of PARAM-2 plots, since 
they represent another three dimensional surface projected into the 
alpha beta parameter plane. In that case, however, each omega curve in 
the PARAM-5 graph represents a different PARAM-2 plot. Consider for a 
moment the W=20 rad/sec curve of figure 4-23. This curve depicts the 
topological relief of a cross section of the PARAM-2 magnitude contour 
graph when that cross section occurs at beta = 0.025. Similarly, dif- 
ferent omega curves of figure 4-23 represent such a cross section in 
other PARAM-2 plots. 

The correlation between PARAM-5 and PARAM-1 as well as PARAM-2 pro- 
vides a rather important conclusion. One normally wants to get a good 
indication of system behavior by simultaneously varying alpha and beta. 
It is not necessary to construct a large number of PARAM-1 plots to do 
this. Neither is it necessary to construct a large number of PARAM-2 
plots. By constructing only a few PARAM-1 graphs, as well as two or 
at most three PARAM-2 graphs, one can gain valuable insight into the 
system behavior in the voids by studying the PARAM-5 plots. Thus for an 


initial system design study, one might be inclined to use the following 
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six graphs: two PARAM-1 plots, two PARAM-2 plots, and two PARAM-5 

plots. The constant magnitude values of the PARAM-1 plots would most 
likely be chosen from the bandwidth and resonant peak criteria. The 
constant omega values for the two FARAM-2 plots would normally be dictat- 
ed by the cut-off frequencies. Alternatly, the peak resonant frequency 
might be used as a governing choice. By choosing an intermediate value 
of alpha and beta, the PARAM-5 curves would indicate how other PARAM-1 
and PARAM-2 graphs would have appeared. For illustrative purposes, 
suppose a problem called for a -3 decibel bandwidth of 60 rad/sec, with a 
peak resonance magnitude of no more than 6 decibels around 40 rad/sec. 
Construction of figures 4-9, 4-12, and 4-16, one would be able to con- 
clude the following information: 

1) The tangent locus of the omega curves in figure 4-9 define 
the limiting values of alpha and beta as dictated by the 
resonance magnitude constraint. 

2) The intersection of the tangent locus and the constant 
W = 40 rad/sec curve will determine which values of 
alpha and beta will give this desired resonant frequency. 

3) Looking at figure 4-12, one sees that the alpha and beta 
values determined above will not provide a -3 decibel band- 
width of 60 rad/sec. Instead, it can be seen by interpola- 
tion that the bandwidth will be approximately 57 rad/sec. 

4) Figure 4-16 shows the magnitude relationship at the de- 
sired resonant frequency. Thus by studying this graph, 
one sees that a resonant magnitude of approximately 7 
decibels will resuit from a 60 rad/sec bandwidth. 

Tne design example above has shown 4 common problem which exists for 


most engineering studies. The desired performance measures cannot be 
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attained; instead some compromise must be made. This conclusion was 
derived by simply observing three frequency response parameter plane 
graphs. The advantage of PARAM-5 graphs lies in the ability to corre- 
late performance sensitivities with parameter variations. Even though 
the curves shown in figure 4-23 do not coincide with the beta most likely 
to be selected, they do illustrate the behavior of magnitude and omega 

at a neighboring beta value. Thus it is possible to project that in- 
formation onto one of the PARAM-1, PARAM-2 graphs and extrapolate the 
data thus obtained. Additionally, one can do the same with the PARAM-5 
curves shown in figure 4-24. This figure shows the magnitude and omega 
variations as a function of beta for a constant alpha also near the 

point of interest. In effect the curves in figure 4-24 reflect the be- 
havior about the peak resonant point of W=35 rad/sec. It is seen that 

a variation in beta will increase the resonant magnitude for smaller beta 
values, with a corresponding decrease in magnitude for higher beta values. 
Figure 4-24 shows another important factor. At beta = 0.036 one sees that 
a flat frequency response of 0 decibels will be obtained for a bandwidth 
of approximately 35 rad/sec. This can be verified in figures 4-11 and 
4-25. 

Figures 4-25 and 4-26 have been included to show more vividly how the 
interpolation of PARAM-5 plots might be made. Consider the PARAM-5 graphs 
as topological portrayals of omega. Figure 4-25 shows how these omega 
contours change from those of Figure 4-23 as beta increases. Note that 
all omega curves are reduced in magnitude when beta is increased. Note 
further that a reduction in the number of omega curves will result in a 
loss of important information. The eight curves shown in figure 4-25 
might lead the observer to conclude that magnitude decreases continually 


as alpha is decreased. Only by careful inspection of the Ww = 40 rad/sec 
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curve can one see that this is not the case. By inspecting the complete 
array of curves in figure 4-23, however, one sees a resonance trend for 
negative alpha values. In this figure, resonance occurs at approximately 
55 rad/sec. Using this information, another PARAM-2 graph can be visual- 
ized for constant W = 55 rad/sec. The curves would be similar in shape 
to those shown in figure 4-15 and 4-16, but smaller in size. The center 
of these elliptical curves would occur at beta = 0.036 and alpha = -0.003. 
The analogous procedure could also be employed using either figures 4-24 


or 4-26. 


PARAM-6 

The discussions of PARAM-3 through PARAM-5 have shown how various 
graphical configurations result from taking cross sectional cuts in 
the PARAM-1 and PARAM-2 frequency response parameter plane curves. The 
curves shown in PARAM-6 depict another such cross sectional projection. 
Consider first the curves shown in figure 4-27. They correspond to 
constant beta curves mapped onto the alpha-omega plane for a given pre- 
selected magnitude value of 12 decibels. Accordingly, this graph corres- 
ponds to the PARAM-1 plot shown in figure 4-8. The full utilization of 
the entire set of graphs will now become clear. Effectively the PARAM-6 
graphs show the third variable of the curved elliptical surface which 
was projected onto the parameter plane in PARAM-1. In this manner, the 
third variable, omega, now becomes one of the axis variables while cross 
sectional cuts of the three dimensional alpha-beta-omega surface are ob- 
tained by slicing down along constant beta lines or constant alpha lines. 
The latter cross sections are pictured in figure 4-28. Thus when the 
corresponding PARAM-1 and PARAM-6 graphs are viewed simultaneously, one 


can visualize the three dimensional alpha-beta-omega surface in its true 


15 


form. The illustrative example of an elliptically shape bull's horn is 
applicable in all but one facet which will now be corrected. It was 
originally stipulated that the omega axis passed down the conical surface, 
and was rotated and shifted slightly as magnitude decreased. This stipula- 
tion was made to simplify the visualization of the surface. Inspection of 
the PARAM-6 curves will show a slight modification of this illustration. 
The omega axis does not pass down the surface, but instead projects per- 
pendicularly through the origin of the PARAM-1 plots. In this manner it 
is seen that the surface has shifted slightly, and been rotated in the 
three dimensional alpha-beta-omega space! Now the illustration is cor- 
rect in its entirety. Analyzing the corresponding PARAM-1l and PARAM-6 
graphs gives a complete description of the portrayed surface in three 
dimensions. Inspection of any one of the PARAM-2 plots reveals the 

change in shape of this three dimensional surface as the fourth dimen- 
sional parameter, magnitude, is varied. Figures 4-29 through 4-32 are 


included for amplification and verification of the above results. 
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GLAVIS. ~PARAM 4. 3RD ORDER I-A F-B., MAG CURVES 


GRDINATE=A_PHA, ABCISSA=OMEGA, CONBETA=0.036 


FIGURE 4-21 
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GLAVIS, PARAM 5, 3RD ORDER TVA F.B.,OMEGR CURVES 
GRDINATE=ALPHA, ABCISSA=MAG ., CONBETA=0. 025 


FIGURE 4-23 
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FIGURE 4-24 
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GRDINATE=ALPHA, ABCISSA=MAG . CONBETA=0.036 


FIGURE 4-25 
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GLAVIS, PARAM 6. 3RD GRDER [~A F-B--, BETA CURVES 
SRDINATE=ALPHA. ABCISSA=OMEGR. CONMAGNI TUDE=12DB 
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ORDINATE=BETA . ABCISSA=OMEGA. CONMAGNI] TUDE=12DB 


FIGURE 4-28 
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GLAUIS. PARAM 6, 3RD, ORDER 178 7-6. see eRUES 
ORDINATE=ALPHA. ABCISSA=OMEGA. CONMAGNITUDE=6DB 


FIGURE 4-29 
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FIGURE 4-30 
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GLQUIS. PARAM 6. 3RD GRDER T/A F.B.. BETA CURUES 
JRDINATE=ALPHA, ABCISSA=GMEGA, CONMAGNITUDE=3DB 


FIGURE 4-31 
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FIGURE 4-32 
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Sk CONCLUSIONS 

Section 4 presented an exhaustive study of the application of fre- 
quency response parameter plane curves to one particular system. -The 
selection of the third order system was based on the ease of interpre- 
tation of each frequency response graph discussed. Extensions to systems 
of higher order will undoubtedly change the graphical configurations to 
some extent. The degree of this variation was smaller than expected for 
a similar fourth order system. Figure 5-1 shows the particular fourth 
order tachometer and acceleration feedback system which was studied. 

Figure 5-2 presents the system stability performance curves in the alpha- 
beta parameter plane. Figures 5-3 through 5-12 depict the associated 
frequency response curves. Comparison with the third order system discuss- 
ed in section 4 shows strikingly similar results. 

Several concluding comments are noteworthy concerning the use PARAM-1 
through PARAM-6. The alpha - beta parameter plane provides more informa- 
tion than any other representation for system studies involving the simul- 
taneous variation of two parameters. Several different PARAM-1 plots 
are considered to be most useful for initial system studies, when they 
are accompanied by their corresponding PARAM-6 graphs. In this manner, 

a three dimensional space of alpha, beta, avd omega can be visualized for 
each of the constant magnitudes used. The additional incorporation of one 
or two PARAM-2 graphs will describe the effect of a variation in the fourth 
dimension. Stated in other words, the information in PARAM-2 shows how the 
three dimensional PARAM-1, PARAM-6 surface varies in magnitude. PARAM-3 
and PARAM-5 curves yield a more sensitive portrayal of the magnitude varia- 
tions; however these are restricted to a single parameter variation. They 
are useful for accurately presenting a maximum peak resonance magnitude 
when this information is not available from the peak resonant tangent 
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locus of PARAM-1. PARAM-4 effectively shows a series of Bode curves 

as a function of one system parameter. As such they would be advantag- 
eous for intermediate studies, when the final adjustments are being 
contemplated. Finally, the familiar Bode plot representations of PARAM-7 
show the most accurate frequency response, and are therefore most useful 


in the terminal stages of design. 
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FIGURE 5-4 
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GLAVIS., PARAM 6 47H ORDER 17 F«Bew AIPA CURYES 
ORDINATE=RALPHA,. ABCISSA=OMEGA, CONMAGNI TUDE=6DB 


FIGURE 5-12 
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GO TO 448 


436 WRITE(696) BETAsCONST,ZMAG y,OMEGA 
GO TO 448 
437 WRITE(6¢6) BETA,CONHET »-ZMAGSQ(LOM) »-OMEGA 
GO TO 448 
438 wRITE(6,6) CONALP, 8ETA,ZMAGSQ( LOM) - OMEGA 
GO TO 448 
439 wRITt(6¢6) BETA, CONBET »ZMAG,OMEGA 
GO TO 448 
440 WRITE( 6,6) CONALP,BETA,ZMAG,OMEGA 
GO TO 448 
441 wRITE(5,6) BE TA,CONST ¢CCNMAG »OMEGA 
O TO 448 
442 conn ae CONST »8ETAsCONMAG » CGMEGA 
443 coe ie ZALPHA(LCG™) -CONBET »YMAGDB, OMEGA,TANG 
444 Gite cna. CONALP,28ETA(LOM) »Y4AGDBe OMEGA, TANG 
ile 
oe Wee eaee +6) CONST ,BETA,ZMAG,OMEGA 
44% 
eo aa ARRAY STORAGE SECTION€ €e& eee eae KEE KE RE EKER 
45 =[P+] 
Y(IP)= PRET 
452 GO TO (453 534460 e470 946074709475 41000) »NPARAM 
453 X(1P)= CON 
GO TO $00 
450 X(1P)=ZMAG 
GO TO SCO 
470 X(IP)=\JMEGA 
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J= 
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La 
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usr 
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ee. 


LOGIC OUTLINE OF PROGRAM PARAMS 
Read number of runs submitted. 
Set run number equal to l. 
Read input data for appropriate run. 
Print input data for appropriate run, if requested. 
Calculate graph scales and location of origin. 
Set NPARAM equal to the first PARAM value. 
Select initial values for appropriate PARAM. 
Print graphical information for appropriate PARAM, if requested. 
Construct logarithmic grid if required. 
Set curve number equal to l. 
Select curve value for appropriate curve number. 
Set abcissa variable to minimum allowable value. 
Set number of graph points equal to 0. 
Compute numerator and denominator polynomial coefficients. 
Check for possible division by 0. 
Calculate ordinate value, using the quadratic formula for PARAM-1 
through PARAM-6. 
Check if calculated ordinate value is within graph limits. If yes, 
continue; if no, go to statement 21. 
Increase number of points by l. 
Print computed values. 
Store values necessary plotting. 
Increase abcissa variable by pre-computed increment. 
Check if abcissa variable exceeds maximum graph limit. If yes, 


continue; if no, go to statement 14. 
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23. 


24. 


25. 


26. 


a7. 


26. 


29. 


30. 


Se. 


Check if any points were found within graph limits. If none, 

so state and continue without plotting. If one, so state 

and continue without plotting. If more than one, continue and 

plot the graph points. 

Check if both positive and negative coefficients were used in the 
quadratic formula solution where applicable. If yes, continue; 

if no, change sign and return to statement 12. 

Increase curve number by 1. 

Check if curve number exceeds the maximum number of curves stipulated 
in input data. If yes, continue; if no, return to statement ll. 
Draw graph and plot computed curve data on same. 

Increase parameter number by l. 

Check if parameter number exceeds maximum. If yes, continue; if no, 
return to statement 7. 

Increase run number by 1. 

Check if run number exceeds maximum. If no, return to statement 3; 


if yes, stop. 
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FLOW CHART FOR PARAMETER SIX 


Read Input Data 


IWRITE 


OoT+ 


Print Input Data 


Calculate X+Y Scales 


IWRITE 
——— 


(Print Graphical Information 


Magnitude ee Magnitude 


Squared-—~<——< IMAG-1 2 Actual 
Actual f | 


Magnitude} in Decibels 









4 Constant 
Cee a 
Curves | 


Numcurves + NALPHA 


DO ICURVE = 1, NUMCURVES 


TRIG — > —>—- 


au 
Gal >= 
+ 





CURVE = ZALPHA(ICURVE) 
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